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16. Abstract 


Significant densities of doubly charged ions exist in the discharge chambers of electron bom- 
bardment ion thrusters. These ions are undesirable because they are a major plasma constituant 
effecting the sputtering damage which limits thruster lifetime. It would be desirable to reduce their 
density while maintaining good thruster performance. The development of a model which predicts the 
doubly charged ion density is discussed. The accuracy of the model is shown to be good for two dif- 
ferent thruster sizes and a total of 11 different cases. The model indicates that in most cases more 
than 80% of the doubly charged ions are produced from singly charged ions. This result can be used to 
develop a much simpler model which, along with correlations of the average plasma properties, can be 
used to determine the doubly charged ion density in ion thrusters with acceptable accuracy. Two dif- 
ferent techniques which can be used to reduce the doubly charged ion density, while maintaining good 
thruster operation, are identified as a result of an examination of the simple model. First, the 
electron density can be reduced and the thruster size then increased to maintain the same propellant 
utilization. Second, at a fixed thruster size, the plasma density, temperature and energy can be re- 
duced and then to maintain a constant propellant utilization the open area of the grids to neutral 
propellant loss can be reduced through the use of a small hole accelerator grid. The reduction in 
the values of the plasma properties causes a decrease in the doubly charged ion density. 
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INTRODUCTION 

Electron bombardment ion thrusters are presently being considered 
for use in deep space probes and for satellite stationkeeping functions. 
These devices which have the advantage of very high specific impulses 
also have the attendant disadvantage of low thrust densities. This low 
thrust characteristic necessitates thruster operation for long periods 
of time in order to accomplish typical missions (of the order of 10,000 
hours). Thruster lifetimes can be limited by the erosion of ion 
chamber component parts with most of this erosion (or sputtering dam- 
age) being caused by doubly charged ions. Long thruster lifetimes 
therefore require control of doubly charged ion densities. 

Many experiments could be performed to determine how the thruster 
should be operated so that thruster performance would be good and the 
doubly charged ion density would be reduced to an acceptable level. 
However, these experiments would have the disadvantage of being time 
consuming and costly and the results might be applicable to one size 
and type of thruster only. Instead a theoretical model could be de- 
veloped which would accurately predict the doubly charged ion density 
over a wide range of conditions and thruster sizes. This model could 
be applied at low cost to determine the factors affecting the doubly 
charged ion density and how they should be adjusted to reduce the 
double ion density. It should also indicate what effects these changes 
would have on thruster performance. Such a model has been developed 
for electron bombardment ion thrusters and has been verified experi- 
mentally for thrusters which use mercury propellant. A discussion of 
the model's development and verification is presented in this paper 
along with some results and conclusions based upon the model. 


THRUSTER OPERATION 


Many of the assumptions and approximations used in the development 
of the model are based upon a knowledge of thruster operation. This 
section will briefly discuss thruster operation so that the develop- 
ment of the model 1 r. the next section will be more easily understood. 

An ion thruster has two basic tasks to perform: 

1) Ionization of the neutral propellant atoms. 

2) Acceleration of the ions to high velocities producing thrust. 
These two topics will form the basis for the discussion of thruster 
operation. 

Figure 1 shows a schematic for a typical electron bombardment ion 

thruster. The specific type shown has a strongly divergent magnetic 

field which is presently the most common type. However, the operation 

. .. ( 1 ) 

of all types of electron bombardment ion thrusters is very similar 
Electrons are emitted from the cathode and are drawn toward the anode 
which is biased 30-40V positive with respect to the cathode. These 
electrons (called primary electrons) are injected into the primary 
electron region with an energy slightly less than that associated with 
the 30-40V anode voltage. Electrons in this region are kept from going 
irunedi ately to the anode by a magnetic field set up between the cathode 
pole piece and the anode pole piece but collisions eventually facilitate 
electron diffusion across these magnetic field lines so that they can 
be collected by the anode. As a result of the magnetic field contain- 
ment the electron density is much higher (10' 1 cm - ) within this region 
than it is outside of it. The primary electron region's boundary is 
defined by the surface of revolution of the critical (magnetic) field 
line and the screen grid. Because the strength of the magnetic Meld 


MAGNET 
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is fairly low within it, a fairly uniform electron density exists 
throughout the entire region. 

Neutral propellant (e.g. mercury) is injected into the upstream 
end of the discharge chamber at low pressure (10~ 4 torr). Most of the 
interactions between electrons and neutral propellant atoms take place 
in the primary electron region because higher electron densities and 
energies exist there. Electrons bombard the neutral atoms occasionally 
knocking an outer shell electron loose from the atom forming a singly 
charged ion. The production of a mercury single ion requires more than 
10 eV of energy from the incident electron. This electron and the 
ejected electron then share the remainder of the energy which the in- 
cident electron carried originally. This reaction results in the re- 
placement of one high energy electron with two lower energy electrons 
which rapidly randomize with similarly generated electrons to form a 
Maxwellian electron group. Ions are extracted from the plasma through 
holes in the screen and accelerator grids as a result of the large po- 
tential difference applied across these two grids. The rate of ion 
loss through the grids times the ionic charge is called the beam 
current. 

Electron bombardment of atoms and ions also produces doubly 
charged ions. Many of these ions are extracted from the discharge 
chamber by the grids, however, some of them go to the walls. As ions 
near the walls (the cathode pole piece, screen grid, etc.) they are 
accelerated to high velocities by an electric field that exists at the 
plasma boundary. When these high velocity ions strike the walls they 
can knock atoms loose (sputter atoms) from the walls of the discharge 
chamber. The energy that doubly charged ions possess upon striking the 
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walls is twice that of singly charged ions, therefore the sputtering 
damage caused by a double ion is much greater than that caused by a 
single ion. Double ions are thought to cause most of the sputtering 
damage even though their density is typically an order of magnitude 
less than that of the single ions. 


THEORETICAL MODEL 


1 


Introduction 

In order to develop a simple model for determining the double ion 
density in the discharge chamber only those ionic and atomic species 
which were considered significant in determining the double ion density 
were included. The significant species were selected as those which 
have substantial electron impact cross sections of formation over the 
electron energy range of interest so that large numbers of these ex- 
cited atoms or ions will be produced. These states also have suffi- 
ciently long effective lifetimes so that they can participate in pro- 
duction processes before they decay. Only those reactions which lead 
directly or indirectly to the production of double ions were included. 

Figure 2 is a discharge chamber reaction schematic showing these 
dominant species and the reactions in which each specie can participate. 
The model has been developed for thrusters using mercury propellant but 
the general procedure is valid for thrusters using other propellants. 

The symbols used in Figure 2 represent the following species: 

Hg° - neutral ground state mercury 

Hg m - metastable neutral mercury (6 3 P 0 and 6 3 P 2 states) 

Hg r - resonance state neutral mercury ( 6 3 P i and 6 1 P j states) 

Hg + - singly ionized ground state mercury 

Hg m+ - singly ionized metastable mercury (6'D^ and 6-'D,^ states) 
Hg ++ - doubly ionized ground state mercury 
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The arrows in Figure 2 indicate the various interaction routes 
considered in this analysis. Three different types of reactions are 
indicated in this figure. The first type of reaction occurs when an 
electron interacts with an atom or ion producing a more highly excited 
specie. This reaction is indicated in Figure 2 by an arrow going from 
one specie to another more highly excited specie (e.g. the production 
of double ions from single ions). The production of a highly excited 
specie also represents a loss mechanism for the less excited specie. 

The reverse reaction in which, for example, an ion captures an electron 
is improbable because the reaction requires three bodies to simultane- 
ously col 1 ide. 

The second type of process is that of an atom or ion going to a 
plasma boundary. Such a boundary could be either the discharge 
chamber wall on which the atom or ion would be de-excited or it could 
be a grid aperture in which case the atom or ion would be extracted 
from the discharge region. In either case this represents a loss rate 
for any of the excited states. These losses to the boundary are in- 
dicated in Figure 2 by the dotted lines to the wall of the chamber. 

The large arrow back to the neutral ground state represents the 
resupply of neutral ground state atoms either from the walls or from 
the propellant supply system. 

A third type of reaction shown in Figure 2 is relevant only to the 
two resonance states. The resonance states differ from metastable 
states in that they have a very shot t lifetime before they de-excite 
spontaneously by emit ting a photon of light. However, the energy of 
this photon is -inch that, it is readily absorbed by a nearby neutral 
ground state atom producing another resonance state atom. Since the 
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transport time of the photon is small compared to the excited state 
lifetime the excited state can be considered to exist continuously. 
Eventually the photon can diffuse to a boundary where it will be lost; 
this is equivalent to the loss of a resonance state atom. This loss 
mechanism is represented in Figure 2 by a dotted line conveying a 
photon to the wall and a branching line going from the resonance atom 
to the neutral ground state atom. 

Figure 2 also shows the dominant routes for the production and 
loss of all of the excited atoms and ions considered. For example, 
ground state single ions can be produced as a result of electron 
bombardment of neutral— ground-s.tate, resonance state, and metastable 
state atoms and these single ions can be lost as a result of single 
ion migration to the plasma boundary and the production of metastable 
single ions and double ions by electron bombardment. 

When equilibrium conditions exist in the discharge chamber the 
rate of production of each specie must equal its loss rate. If, for 
example, the production rate of single ions increases, the single ion 
density must also increase to keep the loss rate (which is directly 
proportional to the single ion density) equal to the higher production 
rate. This example illustrates the fact that the equilibrium density 
of any specie is determined by the associated production and loss rates. 
If equations determining the production and loss rates could be de- 
rived, these equations could then be solved for the equilibrium density 
of any specie under consideration. The remainder of this section is 
concerned with deriving equations for the production and loss mech- 
anisms indicated in Figure 2 and then solving these equations for the 
equilibrium densities of the various states. 
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Electron Bombardment Reactions 

The first reaction to be considered is the one which produces ex- 
cited atoms or ions by electron bombardment from less excited atoms 
or Ions. The total rate of production of any specie y from specie a 
(and hence the loss rate of a due to this reaction) is given by: 

r y = I J n a Y ( E) v (E) dn dV- 0) 

01 Plasma Volume E=0 a “ 

where n is the density of specie a at some point r in the plasma, 

a 

o y (E) is the cross section for the production of y from a at the elec- 
tron energy E (2_6) , v g is the electron velocity at energy E, dn g is the 
density of electrons with energies between E and E + dE at r, and dV- 
is the infinitesimal volume element. The distribution of electrons 
over the energy spectrum of an ion thruster was assumed to be composed 
of a Maxwellian electron group which is described by a temperature 
(T - eV) and a density (n my - cm" 3 ) and a monoenergetic group 
(primary electrons) which is described by an energy (fp^ and a 

density (n - cm" 3 ). This type of electron distribution is generally 
P r (i,7) 

accepted as appropriate for electron bombardment thruster plasmas. 

Substituting this electron distribution into Equation (1) and com- 
bining terms to form new functions results in the following equation. 


R y = 




/ n X n 

Volume 


pr 


p : <y> 


+ n 


mx 


(T )] dV- 
' mx /J 
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p Y (■; ) = V (• 
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) 0 Y (-' 



where 


(3) 
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and 


E=« dn (E) 

qX*> ■ l . »I< E ) ».(E) -¥■ — 


E=0 


(4) 


mx 


"^ n mx^^ n mx^" 1S the Maxwellian distribution function and the other 
terms are as defined previously. 

Where possible the cross sections (ojj required for Equations 
(3) and (4) were selected from experimental data. ^’^>4) eX p er -j_ 
mental data were not available, theoretical cross sections were either 
obtained from the literature^ or calculated using the Gryzinski 
approximation.^) The Gryzinski approximation was modified for the 
cases of the metastable single ion production cross sections to re- 
flect the significant value of the cross sections near the threshold. 

The cross sections used are presented in Figure 3 along with references 
indicating their origin. 

Using integral equations like Equation (2) in the model would be 
inconvenient because it would then be very difficult to solve for the 

density of specie a (n ) since n appears within the integral. For 

a a 

this reason it would be desirable to convert Equation (2) into a simple 
algebraic equation. Fortunately the plasma is fairly uniform in the 
primary electron region which is where most of the reactions take place. 
This suggests using average properties in Equation (2) to obtain the 
following result. 



(5) 


The asterisks indicate volume averaged quantities and V is the volume 
of the primary electron region. 
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In order to use Equation (5) to evaluate all of the production 

* * * * 

rates the volume averaged plasma properties (n pp , n mx , c pr . T^) must 
be determined. Comparing Equations (2) and (5) the following defini- 
tions of the volume averaged properties must apply. 


★ 

n 

a 


* 

n 

a 


pr 


mx 


P I<V> * * 

WO * ■ 


/ n n 
Volume a 

I n 
Volume a 


P Y U ) dV- 
pr a v ^pr' 

n mv Q Y (T mv ) dV- 
mx mx' 


( 6 ) 

(7) 


These two equations show that the volume averaged plasma proper- 
ties will be weighted in some manner. In order to evaluate the 
integrals, species a and y must be chosen. The only specie density 
(n^) that can be determined readily is the single ion density because 
plasma neutrality requires it to be approximately equal to the electron 
density. Specie y must also be chosen in order to determine what P Y and 
Q y to use. Figure 2 shows only two choices are possible— the singly 
ionized metastable states and the doubly ionized ground state. Since 
the whole purpose of the model is to determine the double ion density, 
specie y was chosen as the doubly ionized ground state. Using these 
choices for species a and y , Equations (6) and (7) were rewritten in 
the following form where the electron density (n g ) has been used to 
approximate the single ion density. 


n e V * * L._."e V P + + <V d * 


■pr 


Volume 6 pr 


pr ' 


( 8 ) 


n e Q+ (T|llx) V ‘ = Volume^ " mx Q+ ^ ** 


(9) 
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The volume averaged values of the primary electron energy (c ) 
and the Maxwellian electron temperature (T ) were defined as shown 

1 1 1 A 

in Equations (10) and (11). These definitions were chosen because they 
give reasonable values for the properties involved (i.e. these volume 
averaqed values can't be greater than the ppak values a which wss 
possible with some of the other definitions). 

P ++ fr* i - Volume e pr ■ 

p + (5 pr ) - 7 (10) 


) ++ (T* i 

4 U mx ; 


J p 

_ Volume 

pr 

' + ' <j pr' 

/ n 

Volume 

n„ 

dV- 

pr 

/ n e 

_ Volume 

V 

«■ 

/ n e 

n m „ 

dV- 

Volume e 

mx 


(ID 


Equations (8) - (11) along with Equation (12), which says the volume 
averaged electron density is the sum of the volume averaged primary 
and Maxwellian electron densities, can be combined to obtain the 
following definitions of the remaining volume averaged plasma 
properties . 


V + V 


( 12 ) 


f, 
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(13) 

(14) 


( 15 ) 
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This concludes the mathematical development for electron bombardment 
reactions . 

The volume averaged plasma properties must be evaluated in order 
to use Equation (5) to calculate the production rates. The plasma 
properties (np r , n mx , ^p r , T mx ) are measured at many points inside the 
discharge chamber by a Langmuir probe. This data is then used to 
evaluate the integrals in Equations (10), (11) and (13)-(15) numerically 
yielding the needed volume averaged plasma properties. 

Migration Losses 

The second type of process to be considered is that of an excited 
atom or ion going to the plasma boundary. The equation for the plasma 
boundary loss rate of a specie a is given by: 

R = / n v dA (16) 

plasma boundary 

where n is the density of specie « at the boundary, v is its average 

a a 

velocity toward the plasma boundary, and dA is the infinitesimal area. 
For neutral particles assumed to have a temperature equal to the dis- 
charge chamber wall temperature and having a mass "m 0 " the average 
velocity toward the boundary (v 0 ) is equal to one-fourth the average 
thermal speed 


v., 


i &n 

4" 7 -m,, 


( 17 ) 


whore k is the Boltzmann constant. 
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For ions this velocity is determined by the Bohm criterion 
and is given by 


(8,9) 


v 


q 



n n „ 

o ♦n Et > 
mx 


(18) 


where q is the ion charge (coul) and m i is the ion mass (kg), and 

"'mx* n pr and n mx are as de fi ned Previously. 

Since the integral equation used to define production rates has 
been reduced to an algebraic equation in terms of average properties 
Equation (16) should also be simplified in this manner. Because the 
migration loss is a surface phenomenon, however, it is necessary to use 
surface averaged densities and velocities based on surface average 
properties to obtain 


R 


£a 


n s v s A 

a a 


(19) 


where A is the total surface area of the primary electron region and 
the superscript "s" designates values based on surface averaged prop- 
erties. Equation (19) could be made more convenient for use in the 
model if it were based on volume averaged properties as Equation (5) 
is. Equation (19) was for this reason rewritten in terms of volume 
averaged densities as follows 


R. = n* v* A/F (20) 

9. a a a a ' ' 

★ ★ * ★ 

where v is v (T , n / n ) and F is a plasma uniformity factor given 
t t iiiA P' m X ; t 

by Equation (21) which relates the volume and surface averaged 


density-velocity product. 
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F 

a 


* rt 

n v 

a a 

n s v s 

a a 


* * 


( n v dA 

J a a 

plasma boundary 


/A 


( 21 ) 


This concludes the mathematical development for the migration 
loss of excited species, but some additional discussion and quantifi- 
cation of the uniformity factor F is necessary before it can be used 
in Equation (20). The evaluation of for neutral excited states is 
difficult because it is difficult to measure their densities. The 
migration of excited neutrals to the plasma boundary is however a 
minor loss mechanism compared to the losses due to the conversion of 
neutral excited atoms into single ions. Therefore F^ for neutral ex- 
cited states can be set equal to unity without introducing a signifi- 
cant error into the total loss rate calculation. In the case of ions, 
however, migration to the boundary is a major loss mechanism and F^ 
must be evaluated in order to obtain accurate results. For the singly 
ionized ground state the approximation, n + = n g , can again be used in 
order to evaluate F + using Equation (21). The uniformity factor for 
the singly ionized metastable states was set to unity since these states 
have a very minor effect on the double ion density. The determination 
of F ++ is based on the observation^ 0 ^ that the volume averaged double 
ion density (n* + ) is proportional to the volume averaged electron 
density squared (n*) 2 . It has been assumed that this proportionality 
holds locally and this results in the following definition of F ++ . 


/ n n v ++ . dA j /A 
plasma boundary 


( 22 ) 
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Photon Diffusion Losses 

The third type of process to be considered is the loss of res- 
onance state atoms due to photon diffusion to the walls of the dis- 
charge chamber. From diffusion theory the rate of photon loss across 
any plasma boundary and hence the rate of resonance state atom loss by 
this mechanism is given by the equation 

R £r = / DAn p dA = DA[An p ] . (23) 

plasma boundary 

In this equation n p is the photon density and D is the photon diffusion 
coefficient which is given by 


D = 


mv c > 2 


(24) 


V in this equation is the average lifetime of the resonance state 

atom, n* is the neutral ground state atom density, and is the cross 

0 (11 

section for absorption of the photons by neutral ground state atoms. 

The second equality in Equation (23) reflects the fact that average 
properties are being used in this analysis. The photon density has 
been assumed constant up to a point one photon mean free path from the 
boundary. From this point the density is assumed to decay linearly to 
zero at the boundary. This assumption yields the following conserva- 
tive estimate for the photon loss rate 

R = DA n* hr ( 25 

? r p t 

where < f is the mean free path for photon absorption (-Jr") • This 

1 Vc 
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approximation is valid when the photon mean free path is much less than 
the characteristic dimension of the plasma, a condition that is readily 
satisfied for this case where the photon mean free paths are very small 
U* - .1 cm). 

Since the neutral density is assumed uniform over the discharge 
region the photon density profile is similar to the resonance state 
atom density profile and the following approximation between the photon 
and resonance state atom density at any location in the plasma applies: 


n p = n r !•; . (26) 

P is a proportional i ty constant that can be thought of as the ratio of 
the probability that a photon will be "free" in the plasma to the prob- 
ability that it will be "bound" forming a resonance state atom. This 
ratio of probabilities can also be expressed as the ratio of the 

average lifetime of a free photon ( — * — ) t0 the resonance state atom 

cn a 
o c 

lifetime (i). Therefore, n is given by: 


1 



(27) 


where c is the speed of light and the other quantities have already 
been defined. Combining Equations (25), (26 and { 21 ) one obtains the 
following equation for the loss rate of resonance state atoms due to 


pho ton d i f fits i on : 
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Determination of Specie Densities 

The equations derived so far in this section can now be combined 
to determine the equilibrium density of each specie included In the 
model. Equations of the form of (5) and (20) — and (28) for the case 
of resonance state atoms -- along with the values for the volume 
averaged plasma properties and the plasma uniformity factors can be 
used to determine the rates of production and loss for each specie in 
the plasma. The steady state density of these species can then be cal- 
culated by equating their total production rates to their total loss 
rates. For example, the 6 3 Pq metastable atom density is determined by 
equating the production rate of this metastable atomic state from 
neutral ground state atoms to the sum of the associated loss rates due 
to 1) migration to the wall, 2) production of single ions, and 3) the 
production of double ions, that is 


★ ★ 

n Tn \ j. n m rr* \i u _ n m V m ” 

o L pr V mx Q o (T mx )J * ' TT 


m 


(29) 

+ n m^ n pr P m^pr + n mx V^mx^ ^ + n m ^ n pr P m ^pr^ + n mx ^m ^ 


mx' 


Solving this for the metastable atom density ratio one obtains 

* k 

7" = P o^’pr ^ + n mx Q o (T mx^ 1 {v/aV + ^pr p in^V^ 
o m 

+ n mx V T mx^ + ^- n pr P m ^pr- + n mx ^ T mx^} 


(30) 


★ + 
where n m is the volume average metastable state density and v m is the 

average velocity of metastable neutral atoms toward the boundary. A 

similar type of equation can be derived for each of the other excited 

states but. they are all as complex or more complex than Equation (30). 
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For example, the equation for the doubly ionized ground state density 
has eight terms in the numerator. Each of these terms has the same 
form as the bracketted quantities in the numerator of Equation (30). 

The quantity V-/A in the denominator of Equation (30) has an in- 
teresting physical interpretation. It is contained in a term which 
represents the loss rate per unit volume of metastable state atoms to 
the plasma boundary. This term shows the manner in which the size and 
shape of the primary electron region enters into the model. For a 

large thruster V-/A will be large and an ion or excited neutral must, on 

the average, travel great distances to reach the plasma boundary, and 
so the loss rate per unit volume of these species will be small. For a 
small thruster ¥•/ A will be small and on the average the ions and ex- 
cited species are near the boundary and can reach it readily resulting 

in a large loss rate of these species per unit volume. 

At this point only the relative density of each excited specie 
★ * 

(n /n ) can be calculated. However one additional fact can be added to 

★ 

the model, the requirement that the plasma be neutral (i.e., n g 

k k k 

= n + + n m+ + 2n ++ ). This requirement when added to the relative density 
equations of the ionized states implies unique single, metastable single 
and double ion densities. These in turn imply a unique neutral atom 
density and hence a unique density for each specie considered in the 
analysis. One must however iterate to arrive at these densities be- 
cause a neutral ground state atom density must be assumed initially to 
determine photon loss rates from Equation (28). At the conclusion of 
the analysis then the calculated ground state atom density must agree 
with the assumed value. 
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A computer program has been written which calculates the densities 
of all the species considered in the model. The densities are calcu- 
lated by using relative density equations similar to Equation (30) and 
the plasma neutrality condition. The input needed to make these calcu- 
lations includes the volume averaged plasma properties, the plasma 
uniformity factors and the volume-to-surface area ratio of the primary 
electron region. A listing of this computer program entitled "HG" is 
included as Appendix A. 


EXPERIMENTAL PROCEDURES AND RESULTS 


The model developed in the previous section will predict the 
specie densities if the volume averaged plasma properties, uniformity 
factors and geometric quantities, which are collectively called the 
input parameters, are known. In order to verify the accuracy of the 
model, data must be gathered so that the model's input parameters can 
oe determired. These input parameters can be used by the model to pre- 
dict the specie densities which can then be compared to the measured 
densities to determine the model's accuracy. The model's accuracy has 
been determined by comparing the measured and predicted double ion 
densities since the double ion density is the model's main concern. 

In order to test the accuracy of the model over a wide range of 
conditions data were used from different thrusters and operating con- 
ditions. The 15 cm diameter SERT II thruster was operated with two 
different grid sets and at three different power levels in each of 
these configurations. Data were collected at each condition allowing 
the accuracy of the model to be verified at six different points. Data 

for the 30 cm diameter thruster were also obtained from Hughes Research 
( 12 ) 

Laboratories' ' so that the model could be verified over a wider range 
of thruster sizes, configurations , and operating conditions. Both 
thrusters have strongly divergent magnetic fields. Their general con- 
figuration and manner of operation have been described in the "Thruster 
Operation" section. More detailed thruster specifications, etc. are 
available in the literature.^’ ^ 

The data gathering procedure for the 15 o:.. thruster will bo used 
to illustrate the general manner in which the needed data were obta ined . 
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Before the gathering of data could begin thruster operation and flow 
rates were kept stable for approximately thirty minutes. This insured 
thruster conditions would change little in the twenty minute period 
during which the data were obtained. Table I lists the conditions and 
configurations at which the 15 cm thruster was operated at C.S.U. along 
with those for the 30 cm thruster as obtained from Hughes Research Lab- 
oratories. This table indicates the changes in configuration for both 
thrusters resulted from using different grid types. The SERT II grids, 
listed in Table I, are flat grids with hole diameters of =.4 cm and in 
operation are separated by a gap of .23 cm. The high perveance dished 
grids are dished slightly to prevent the grids from shorting during 
operation due to their thermal expansion. Their hole diameter is 
smaller (.25 cm) as is their separation gap (.079). More detailed 
specifications for the two grid types can be found in Reference 9. The 
EM (Engineering Model) grids are similar to the high perveance dished 
grids described above. The SHAG (Small Hole Accelerator Grids) grids 
have an accelerator hole diameter that is ^70% of the EM grids' accel- 
erator hole diameter. This smaller hole size reduces the loss of 
neutral propellant. These two grid types are described in more detail 
in References 15 and 16. Table I shows, for example, that the 15 cm 
thruster with the SERT II grids was operated at one condition where the 
amount of current collected at the anode (I arc ) was 1.7A, while the 

voltage difference between the anode and cathode (V ) was 37.2 V and 

ci r c 

the ion current through the grids (I beam ) was 0.258 A. 

The values of the volume averaged plasma properties and the uni- 
formity factors must be known in order to calculate the theoretical 
double ion density. In order to determine the values of these average 
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Table I 

Thruster Sizes, Configurations and Conditions 


Thruster 

Diameter 

(cm) 

Grid Type 

Anode 

Current 

< ! arc-- A > 

Anode 

Voltage 

<v c -- v > 

Ion Beam 
Current 

^beanf 

Mass Flow 
Rate 
(A) 

15. 

SERT II 

1 .00 

33. 

.180 

.310 

15. 

SERT II 

1.70 

37.2 

.258 

.307 

15. 

SERT II 

2.05 

42.6 

.272 

.308 

15. 

Dished 

3.02 

32.2 

.439 

.735 

15. 

Dished 

4.06 

37.5 

.654 

.725 

15. 

Dished 

4.13 

40.4 

.622 

.650 

30. 

EM 

5.0 

37. 

1 .0 

1 .25 

30. 

EM 

7.5 

37. 

1.5 

1.76 

30. 

EM 

10.0 

37. 

2.0 

2.29 

30. 

SHAG 

9.5 

30. 

1.5 

1.74 

30. 

SHAG 

11.7 

30. 

2.0 

2.30 
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plasma properties, the primary and Maxwellian electron densities and 

energies must be determined everywhere in the discharge chamber. The 

plasma properties at some point in the plasma can be measured using a 

Langmuir probe and analyzed using the procedure described in 

Reference 17. The plasma properties at sixteen different points in 

the plasma were measured, for each 15 cm thruster condition listed in 

Table I, using the movable Langmuir probe and associated circuitry 

described in Reference 18. The results of a typical survey (15 cm 

thruster - SERT II grids -37 V anode voltage) are plotted in Figure 4. 

This figure shows the spatial variation of the Maxwellian electron 

temperature, primary electron energy and the primary and Maxwellian 

electron densities in the discharge chamber. The Maxwellian electron 

temperature is seen to average approximately 9 eV over the primary 

electron region defined by the critical field line while the primary 

electron energy averages about 30 eV. The average Maxwellian electron 

density is about 10 11 cm" - while the average primary electron density 

is approximately 10 ! °cnr • over the same region. The electron densities 

and energies are seen to be fairly uniform in the primary electron 

region but drop off rapidly outside this region. Using data similar to 

that plotted in Figure 4, Equations (10)-(15) and (21) were evaluated 

numerically by the computer program "PROP" (listed in Appendix B) to 

obtain the volume averaged properties and uniformity factors for each 

case. The results are listed in Table II along with the volume-to- 

surface area ratio of the primary electron region (V/A) and the 

thruster operating specifications which are reproduced from Table I. 

The average values which resulted from an examination of Tigure 4 

(T 9 eV, 30 eV , n 10 ; 'em” , n 10 cm" ) are seen to 

' mx pr ’ nx nr 










Table II 
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agree well with the volume averaged values listed in Table II 

«». V * 29 - 6 «’ V * 8 ' 4 * 10l0cln_ '* V * 7 '° * 

1 0 9 crn“ 3 ) . F + and F ++ , which are defined as the ratio of the volume 
average ion flux to the average flux at the surface of the primary 
electron region, are seen to have values of 2.1 and 2.5 respectively 
for the case being discussed. The average plasma properties listed 
in Table II are observed to cover a wide range in plasma conditions; 
a situation which is desirable for verification of the model. 

The double ion density inside the discharge chamber must also be 
determined. This can be accomplished indirectly by determining the 
double-to-si ngle ion density ratio (n ++ /n + ) in the discharge chamber 
and the single ion density. The single ion density can be determined 
witn sufficient accuracy by equating the single ion density to the 
electron density. The value of the double-to-single ion density ratio 
can be determined from measurements of the ratio of the double ion 

■f -f 

current to the single ion current in the exhaust beam (I /I), and 


the equation 






J + 2/2" 


( 31 ) 


The quantity 2/2 accounts for charge and Rohm criterion velocity dif- 
ferences between double and single ions. 

The quantity "I ++ /J H " was measured using a mass spectrometer 
The methods used for data acorn's i tion and analysis using such a device 
are described in Reference 20 fen the 15 cm thruster data and in 
Reference 19 for the 30 an thruster. The results obtained are listed 
in the last row of Table II. They .hi.v.fu e.-. ample, a double-to-single 
ien i.urrent. r.i* in o' '7,5 -it ‘hr I ’• ’ • u ‘ n - t grid - 17 V 
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anode voltage case. The general trend observed from these data is 
that an increase in power input (I times V arc ) for a certain 
thruster configuration results in an increase in the ratio I ++ /I + . 
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RESULTS AND DISCUSSION 

The values of the average plasma properties, listed in Table II, 
are observed to vary over large ranges. For example, the average 
Maxwellian electron temperature ranges from a low value of 3.3 eV to a 
high value of 12.2. Similarly the average primary-to-Maxwellian elec- 
tron density ratio varies from 0.02 to 0.50. This large variation is 
considered sufficient to allow a general decision to be made about the 
accuracy of the model. Comparisons of the experimental and theoretical 
values of the double-to-single ion density ratio have been used to 
verify the model's accuracy because this quantity (n ++ /n + ) was de- 
termined experimentally. The theoretical and experimental values of 
the double-to-single ion density ratio are plotted as a function of 
propellant utilization in Figures 5 and 6. The curves labeled 
"THEORETICAL" result from predictions made by the model using the 
"Input Parameters" listed in Table II. The curves labeled "EXPERIMEN- 
TAL" result from measurements of the ratio I /I made using the mass 
spectrometer. The trends exhibited by the THEORETICAL an' EXPERIMENTAL 
curves are very similar and the agreement between the THEORETICAL and 
EXPERIMENTAL values of the double-to-single ion density ratio is good 
for plasma physics work with the average error being 3%. The maximum 
error of 40% is observed at low double-to-single ion density ratios in 
the 30 cm thruster. These error values indicate the model is accurate 
over a wide range of plasma conditions and thruster configurations. 

Since the model has been shown to be accurate in its predictions 
of the double-to-single ion density ratio over a wide range of condi- 
tions there is a distinct possibility that the specie densities and 
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reaction rates used to predict the double-to-single ion density ratio 
are also accurate. The remainder of this section will examine the 
model's predictions of these specie densities and reaction rates. 

These quantities are listed in Table III along with the thruster oper- 
ating variables and the model's input parameters which were reproduced 
from Table II . 

The section in Table III titled "Calculated Normalized Densities" 
lists the model's predictions of the normalized densities of the states 
considered in the model where the normalized density of some specie is 
defined as the specie density divided by the total heavy particle 
density. The sum of the normalized densities for any thruster condi- 
tion should therefore equal unity. Table III shows, for example, that 
the 15 cm diameter thruster operating with SERT II grids at 37 V anode 
voltage would be predicted to have 68a. neutral ground state atoms, 19% 
neutral resonance state atoms, 6.9% singly charged ground state ions 
and .2 % doubly charged ground state ions. The normalized density of 
the single ions agrees fairly well in all cases with the 10% value 
quoted as typical in the literature.^ As expected the neutral 
ground state atoms are the most numerous. 

These normalized density trends can be explained in terms of 
variations of plasma properties. For example, the normalized single 
ion density increases with increasing power input (L . times V , ) to 
the thruster in all cases. This occurs because an increase in the 
values of the volume averaged plasma properties ( auses the ratio ol the 
production rate of single ions to the total neutral density to in. .cease. 
The increase in the ratio indicates a smaller total neutral density is 


\ 
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needed to maintain a specified single ion density and so the normal- 
ized single ion density increases as previously observed. 

The last section in Table III shows the calculated production 
rates for singly and doubly charged ions through the various inter- 
mediate states. These production rates have been normalized by the 
total production rate of the specie indicated. The fraction of the 
associated interactions effected by the primary electrons is indicated 
in parenthesis. For example, at the 15 cm thruster's 37 V, SERT II 
grid operating point, 59% of the single ions are produced as a result 
of electron interaction with neutral ground state atoms and 28% re- 
sulted from electron bombardment of neutral resonance state atoms. 

The neutral ground state-to-single ionic state interactions were in- 
duced by primary electrons 23% of the time and by Maxwellian electrons 
the remainder (77%) of the time. 

Thruster performance is determined primarily by the mechanism for 
the production and loss of single ions. The production of these ions 
is, according to this model, quite dependent on the neutral metastable 
and neutral resonance states which are ignored in most other analyses. 
The manner in which single ions are produced however differs a great 
deal between the two thrusters. In the 15 cm thruster most of the 
single ions are produced as a result of Maxwellian electron bombardment 
while primary electrons are unimportant because of their low densities. 
This indicates that for 15 cm thruster operation the primary electron 
region is the important reactio., region because it is the region where 
high densities of high energy Maxwellian electrons occur. In the 30 cm 
thruster, however, relatively high primary electron densities exist and 
since the Maxwellian electron temperature is low most of the single ion 
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production results from primary electron bombardment. So for the 
30 cm thruster the primary electron region is the important reaction 
region because it contains high densities of high energy primary elec- 
trons . 

Table III indicates in all cases a large percentage of the double 
ions are produced from single ions. This is as one would expect because 
the minimum energy required to produce a double ion from a single ion 
is 18.7 eV while 29 eV is required to produce a double ion from a 
neutral ground state atom. As the power input to the thruster increases 
the number of electrons with energies greater than 29 eV increases 
causing the relative importance of the neutral-to-double transition to 
increase. The least energy is required for the production of double 
ions via singly ionized metastable states, but the densities of these 
states are so low that this production mechanism is unimportant. 
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SIMPLIFIED MODEL 

In the previous section it has been shown that most double ions 
are produced as a result of electron bombardment of single ions. In 
order to simplify the analysis of the "Theoretical Model" section the 
other intermediate states for double ion production can therefore be 
ignored with no significant loss in the accuracy of the double ion 
density calculations. In the simplified model presented here the 
approximation is made that the total rate of production of double ions 
equals the rate of production of double ions from single ions. This 
production rate is given by: 


[ v CO*; c<T mx >]* • 


++, * 


pr 


pr' 


(32) 


The total loss rate of double ions is given by the equation 


"++ 


C V A 
F +t 


(33) 


Equating the loss and production rates and then solving the resultant 
equation for the double ion density one obtains 


. . , cC CO - v C(C>i * 

n ++ " n + 


(34) 


v ++ A 


1 + + 


The approximation n^ = n + can now be used and Equation (18) can be sub- 
stituted for the double ion velocity to obtain the following equation. 
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The double ion density can now be determined for a given thruster 
operating condition using this equation and the plots of P + U pr ) and 
Q^ + (T nix ) found in Figure 7 if the volume averaged plasma properties 
and the uniformity factor F ++ are known. This equation will consis- 
tently predict lower double ion densities than the complete model since 
it ignores the production of double ions from neutral states and the 
singly ionized metastable states, but this error should generally be 
small. The error will be greatest for plasmas with high energy elec- 
trons which can produce double ions directly from neutral states. 

The last section of Table III can be used to determine the magni- 
tude of this error for the 11 cases considered in this study. Since 
the simplified model considers only the single-to-double transition the 
error associated with this approximation can be determined from the 
listed value of the percentage of double ions produced from single ions. 
For example, for the 15 cm-SERT II grid - 37 V anode voltage case the 
percentage of double ions produced from single ions is 78 %. This means 
that the value of the double ion density predicted by the simplified 
model would be 78-' of that predicted by the complete model. Examination 
of Table III indicates the double ion densities calculated using the 
simplified model will agree well with the complete model's predictions 
for all the 30 cm thruster conditions because in these cases the per- 
centage of double ions produced from single ions is greater than 97 ' . 
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In each of these cases few electrons have energies in excess of 29 eV 
(the minimum energy required for the neutral -to-double transition). 

The simplified model will however, according to Table III, yield re- 
sults which are generally low for the 15 cm thruster data (e.g. 30% 

low for the SERT II grid - 42 V anode voltage condition) because in 
these cases sufficiently high Maxwellian electron temperatures exist 
to cause a relatively large percentage of the electrons to have energies 
in excess of 29 eV. 

The most accurate way to determine the values of the average plasma 
properties required in Equation (35) would be to conduct a Langmuir 
probe survey of the discharge chamber under consideration to determine 
the plasma properties at many different points and to then use this in- 
formation in Equations (10) to (15) and (21) to determine average plasma 
properties. The collection of the plasma property data is however 
costly and time consuming. For this reason average plasma property 
correlations were developed. The correlating parameters used are com- 
posed of thruster operating parameters (e.g. I ) and geometric 
properties (e.g. V-/A). Using the Maxwellian electron temperature data 
listed in Table III, for example, one obtains the correlation presented 
in Figure 8. The terms used in the correlating parameter are defined 
in Table III. The correlating parameter used in Figure 8 was deter- 
mined by trial and error. The shape of a curve through the resultant 
data points was picked to match the trends observed in the data points. 
For example, the slope of the curve in the neighborhood of the low 
Maxwellian electron temperature points is seen to decrease. This agrees 
with the trend observed in the data and also agrees with a prediction, 
based on inelastic collisijn cross section data, which says a lower 
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bound on the Maxwellian electron temperature should exist roughly in 
the neighborhood of 5 eV^. 

The correlation for the primary electron energy is shown in 
Figure 9. The correlating parameter contains the quantity n c which 
is the corrected propellant utilization. The corrected utilization 
was used in the correlating parameter, instead of the measured pro- 
pellant utilization, because a better fit of the data points resulted 
from its use. The propellant utilization (n) of an ion thruster 
depends upon the plasma properties, the effective open area for the 
loss of neutral atoms through the grids (A Q ) and the effective open 
area for the loss of ions through the grids (A + ). The propellant util- 
ization is defined by the equation 


n + \ A + 

" + v + \ + "ot v o A o 


* 1 - 


n ot v o A o 
n + v + A + 


(36) 


where n * is the total neutral atom density. The primary energy 
(and other average plasma properties) of a given thruster correlate 
with the propellant utilization as defined above, but correlation be- 
tween grid sets having different values of the ratio A o /A + is poor. 

The problem caused by the utilization’s dependence upon grid sets ran 
be corrected by eliminating the ratio A Q /A + from Equation (36) and then 
substituting in its place the value of the ratio A Q /A + for some stand- 
ard grid set. The resultant quantity is the corrected utilization and 
is defined by the equation 
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The open area for the loss of ions (A + ) from a thruster is proportional 
to the open area fraction of the screen grid ( 4 Equilibrium flow 
theory^ can be used to determine that the open area for the loss of 
neutral atoms (A ) is proportional to the quantity (<j> <j> ) / ( + $ ) 
where is the open area fraction of the accelerator grid. These two 
approximations can be used to define the ratio A Q /A + as follows 






4>. 


<f>. + *. 


(38) 


If Equations (36) - (38) are combined the following result is obtained, 


n c = 1 - .5(1 - n ) 


*s + ^a 


(39) 


where the constant ".5" defines <{>_/(* + ij> ) for the standard grid set. 

Figures 10-12 show correlations for the remaining input parameters. 

These correlations were developed by trial and error in a manner similar 

to that u ? -d to obtain those in Figures 8 and 9. It should be noted 

that the correlation in Figure 11 is for the quantity n [¥-/A]" l,b not 

pr 

★ 

the primary electron density (n ). 

It should be understood that the correlations of Figures 8-12 are 
based on data obtained from strongly divergent magnetic field thrusters. 
The average plasma properties predicted using these figures may be in- 
accurate for other types of thrusters (e.g. multipole or radial field 
thrusters). Therefore Langmuir probe surveys should be made for these 
other types in order to obtain good estimates of the average plasma 
properties and hence accurate predictions of the double ion density. 


9 
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Application of the simplified model can be best demonstrated 
through an example. Consider a 15 cm thruster operating at the condi- 
tions defined by the first section of Table IV. The corrected utili- 
zation (n c ) is first calculated using Equation (39) and a value of 
68% is obtained. Next the correlating parameters are calculated. For 
example, the value of the correlating parameter 


f 1 arc V arc)' 2 „ ^ 

V ‘beam ' arc 




. This value indicates the average 


1 . .533 

used in Figure 8 is 8.3 — ■ j y ! 

cm ' 

Maxwellian electron temperature would be 4.6 eV. The remainder of the 
average plasma properties were determined in a similar manner. The re- 
sults obtained are listed in the second section of Table IV. Using the 
values of the primary electron energy and the Maxwellian electron tem- 
perature one can enter Figure 7 and determine P+ + (22 eV) and Q* + (4.6 eV). 
These quantities, together with the average densities, the uniformity 
factor and the volume-to-surface area ratio for this thruster are then 
substituted into Equation (35) to obtain the double ion density as 
shown in the last section of Table IV. The double ion density calcu- 
lated using the simplified model is 6.2 x 10’ cm"' while the value cal- 
culated using the complete model is 5.2 x 10'’ cm"'. The major reason 
for the discrepancy is that the electron temperature in Table IV 
(4.6 eV) is larger than the value used by the complete model (4.3 eV). 
This higher electron temperature causes Q + (T ) to be too large and 
results in the over-estimate of the double ion density. 
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T 
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Table IV. 

Determination of the Double Ion Density Using the Simplified Model 

Measured Thruster Variables 



(15 

cm Thruster) 

*arc 

= 3. amps 

4> s = .67 

V 

arc 

= 32.2 volt 

♦a " ' 67 

*beam 

= .499 amps 

V-/A = 1.4 cm 

n 

= .68 


n c 

= .68 



Approximate Plasma Properties 


* 

T mx 

=4.6 eV 

* 

V 

= 4.64 

x 10 9 cm -3 

* 

? pr 

= 22. eV 

★ 

n e 

★ 

= 3.51 
* 

x 10 11 cm" 
★ 

F 

= 2.55 

n 

= n - 

n = 3.46 

r ++ 


mx 

e 

pr 


Calculation of the Double Ion Density 
n* = (3.51 x 10 11 cm“ 3 ) 2 (1.4 cm) (2.55) 

TT 

[9.6 x IQ 9 C f — (4.6 eV) (1 .013)]^ 
sec -eV 

x [.013 (.55 x 10" 7 frr) + .987 (.23 x 10' 8 |~)] 


6.2 x 10 9 cm’ 3 


5?. 


An examination of Equation (35) will indicate some general trends 
which should be considered in the design and operation of electron bom- 
bardment thrusters. For example, the double ion density varies linearly 
with the volume-to-surface area ratio. Therefore if two thrusters have 
the same average plasma properties the larger thruster will have a 
higher double ion density. Equation (35) suggests it would be desirable 
to reduce the electron density since the double ion density is propor- 
tional to the square of the electron density. However, making arbitrary 
adjustments in the plasma properties to reduce the double ion density 
may have an adverse effect on other aspects of thruster performance 
which must also be considered. An examination of the effect of elec- 
tron density on propellant utilization will indicate one of the effects 
such an adjustment would have. The propellant utilization previously 
defined in Equation (36), is reproduced below. 


The single ion density (n + ) can be approximated by the electron density 
(„*). The total neutral density (n* t ) is the sum of the densities of 
all the neutral species and can be calculated using the equation 
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where and are 
si ties. The values o 
using equations si mil 
state density can be 


the total metastable and resonance states den- 
f tiie ratios in Equation (41) can be calculated 
( ,r in form to Equation (30). The neutral ground 
calculated using the equation 
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(42) 


k k 

where the ratio (n + /n 0 ) again takes a form similar to that of Equation 
(30). Combining these results into Equation (40) a result of the 
following form is obtained, "f" is a function of the Maxwellian elec- 
tron temperature »„_primary electron energy, primary-to-Maxwellian 
electron density ratio and the uniformity factor F + . The dependence 
of the propellant utilization on the electron density and thruster 
parameters is explicitly shown. 
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One can see that a reduction in the electron density to reduce the 
double ion density will also have the undesirable effect of reducing 
the propellant utilization. However, if some changes in thruster de- 
sign are made along with a reduction in the electron density the 
propellant utilization can be held constant while the double ion den- 
sity is reduced. For example, if a new thruster were being designed 
one might double the volume-to-surface area ratio by making the 
thruster larger than its predecessor. It could then be operated at 
one-half the electron density of the predecessor allowing the propellant 
utilization to remain constant while exhibiting half the double ion 
density in accordance with Equation (35). 

It might also be desirable to reduce the double ion density of a 
certain size thruster while maintaining the same propellant utilization. 
The propellant utilization could be held constant by reducing both the 
ratio A q /A + (which reduces the relative escape rate of neutrals) and 


★ 

the electron density in a manner that keeps the ratio (A 0 /A + )/n e 
constant. According to Equation (35) this would result in a large 
reduction in the double ion density which varies as the square of the 
electron density. The data in Table III for the two 30 cm thruster 
configurations at 1.5 and 2.0 amps beam current can be used to de- 
termine if theory and experiment agree for this method of double ion 
density reduction. The only difference in these two thruster config- 
urations is the open area fraction of the accelerator grid. The EM 
accelerator grid has an open area fraction (<j> ) of 45% while the open 
area fraction for the SHAG accelerator grid is 23%. Both sets have a 
69% open area fraction for the screen grid. The value of the ratio 
A q /A + can be calculated for both grid sets using Equation (38). For 
the EM grids the ratio A Q /A + has a value of .39 while for the SHAG 
grids the value of the ratio is .25. The change from EM grids to SHAG 
grids then allowed operation at a given propellant utilization to occur 
at a lower arc voltage and hence a lower electron densities and energies 
and as a result lower double ion densities. In this particular case 
the double-to-single ion density ratios dropped from 4.4% and 6.0% to 
2.2% and 2.8% respectively at two different utilizations when the SHAG 
grids were used. The theoretical model predicted essentially the same 
quantitative changes. 
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CONCLUSIONS 

A discharge chamber model for an electron bombardment ion thruster 
has been developed which considers metastable, resonance and ground 
state atomic and ionic production and loss mechanisms. The model can 
be used to predict doubly charged ion densities from plasma property 
information. These calculated double ion densities agree with measured 
values to within 40% for low values of the double-to-single ion density 
ratio (n ++ /n + < 2%) and to within 20% for the rest of the data. Cor- 
relations, which relate average plasma properties to thruster operating 
variables such as anode current, can be used to estimate the average 
plasma properties in strongly divergent magentic field thrusters when 
the properties themselves are not available. Singly charged ions are 
produced, according to this analysis, in significant numbers in two 
step processes through intermediate metastable and resonance states in 
addition to direct ionization from the neutral ground state. Doubly 
charged ions are produced predominantly via the singly ionized ground 
state with direct ground state neutral-to-double ion production be- 
coming more significant in plasmas with high Maxwellian electron tem- 
peratures and primary electron energies. A simplified model which 
considers only the singly ionized ground state in double ion production 
can be used to predict double ion densities that agree with the 
complete model's predictions to within 5% when primary electron ener- 
gies and Maxwellian electron temperatures are less than 29 eV and 5 eV 

( 12 ) 

respectively. The recent experimental observation^ that the use of 
small hole accelerator grids in conjunction with lower anode voltages 
provides a means for reducing double ion densities in thrusters, with- 
out degrading performance, is supported by the model. 
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APPENDIX A 


The computer program "HG" , which can be used to predict the den- 
sities of excited atomic and ion states considered in the complete 
model, is listed below. The input parameters needed by this program 
can be approximated using the correlations in Figures 8-12. More 
accurate input parameters can be determined using the computer program 
"PROP", listed in Appendix B, and data obtained from a Langmuir probe 
survey of the discharge chamber. The computer program "HG" uses the 
equations developed in the "Theoretical Model" section and carrys out 
the calculations in the manner suggested at the end of that section. 
Comment cards are included in the listing to indicate what calcula- 
tions, etc. are to be carried in each section. 

Values of the functions P y (f 1 and Q Y (T ) are listed immedi- 

a (Jr ot uiA 

ately after the computer program listing. The particular initial state 
(a) and final state (y) are indicated in the last twenty columns. For 
example, the label "HGM-HG+ 3P0" indicates the initial state for the 
reaction is the 6 3 P Q metastable state and the final state is the singly 
ionized ground state. The first seven cards listed with a particular 
identifying label contain the values for P Y (r ) while the second 
seven list values for Q Y (T mx )« 


ooort 


r nm 


PRO^O&M HG (INPUT. OUTPUT, PUNCH. T APF5»1NPUT , T APF6=0'ITPUT . TAPE8«PUNCHG 
1H.FIIMPU 

COMMON /h/ NP«h M » TF.MP.PftINRG 

COMMON /o/ UTT < 4 , ? c > ) » P|) (* ,25 ) *IK . 

DIMENSION A<8), HOT. CO). 0(3). F(3>. SAOO* ( 2 , 2 , 2 1 ) , 5.K00M < 2,2 , ?M6 
11). ?X01(?,?l>, St-niiA^ll. SK02(2,2l). ^02(2.2), SAOMl (2,2.21), HG 
2 ^EOMl (2.2.21 ) , SX11 m(2»2,<?1). SF 1 1* (2,2,21 ) » Sa12(c,21), 

31). PH (?) , TP(?). SfI u 2(2.2. 21). *#0*0(2). HftlMO (2) . bX (2.21) . S MHG 
42,21) , P2 (P) , TI(P). MS), Mfl), 0(H), H(b). Z(S). H(8), SX°R(2.2.HG 
B21). SPOP(2.?,21) • «XH1(?,?.21)» SF»1 (2.2.21 ) ? FP(2). TM(2), W(2).HU 
6 CO). Tk(2), X Sm ( 2 ) . xSPM(2), XPSM(2). X-ShmiM. FNFF (2) , 

72). TAW (?) . XLAM(2). 0PPI6.6). XNFIO). XL (J.fc.10) . Phfc(6»10>. SP(6MS 
8.6), P i- A ( 7 ) . HA ( 7 ) . PYP ( 30 ) » RN0(4>. SO ( 1 4 . 6 . 1 0 ) . TTltl4.fl), PP<10 *® 
9). PP(10). r XP ( 30 ) . SXl M 2(2»2,2l) , TMP(b). 0T(«). f)rM«) 

DATA xnuTH (1,1) .TTM2.1) .TTI (3.1) .TTI C.l) .TTI ( 5. 1 1 ). TT I ( 16 ,1 ) . TTHG 
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THIS PROG RAM is 5F T UP FOP MERCUPY ONLY 
SET UP nr The P»OfiP*M 

Nn 2 = 1 0 
NNN = 0 

DC 1M J=2*8 

h ( I> = 1 Oh 

101 X(1)«10H 
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PFAO (8*140) (TmP(l) fL-ltNS) 

On TO 10$ 

Ft OENSsTMP(JJ) 

GO TO i 09 
WPTTF (6*139) 

CONTINUE 
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CALCULATION OF METASTA^LE/ NEUTHAL GROUND STATE 


PATIO 


00 U? J=1 *2 
00 110 1=1.2 

DO 110 E-l*21 

<t { I «K) =SXOOM( J. I • K > 

S E < I • K ) *= SEOO M ( J • T • < ) 

CALL YIN TEG (SXtSEt 1 . • SU*4*«HSIJH ) 
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C CAI.CULATF RLSuNANCK/NKH IRAL RATIO 

C 

no 117 J-l . N<J 

IK ( IKLAU.FU. 1) GO TO 116 

no lu 1 = 1.2 
OO IP* <-1.21 

5X ( I .* ) -SXO!M J* I * K 1 
114 SK ( [ ♦ K ) =S6 OK ( J, I ,K) 

CALL. YINTKG (SX«Sb'tl»tXSM(j) «XPSM{J) ) 

On 11^ 1=1.2 
On 11S ksI .21 

SX ( T »k ) = 6X^1 ( j, I .k> 

1 IS SK(I.K) = SFPl ( J. t *K) 

CAU. YlNTF.G (SX *SE» 1 • tXS^M ( j) ,XPSMM(J) ) 

116 CON! T I Ml IF 

nv 0=33 9 63, /XL Am ( j ) 
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SO (3* J J* II ) - ( 9- ( 11 ♦PP ( 2) ) <- 10 0 . 
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x* a- i no ./ ( l (ui , ( i . i.j . it ) ♦m ( 2 . jj . II 1 ) 

SU ( 7 . ,lh II ) It* >x K\ 

XM-^'A«ll,* f .;0‘'('(l) ♦HA 0*0(2) ♦r-Ml 1 ♦HU (21 1 
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^ '/ - 0 v ♦ * '*'(.!) t p s *A ( J ) 


HG 1980 
HG 1990 
HG 2000 
MG 20 10 
HG 2020 
MG 2030 
MG 2040 
MG 20 SO 
HG 2060 
Hr, 20 70 
HG 2060 
HG 2090 
MC, 2100 
MG 2110 
HG 2120 
MG 2130 
HG 2140 
HG 2 l SO 
HG 2160 
HG 2170 
HG 2 IMG 
MG 2190 
HG 2200 
MG 2210 
HG 2220 
HG 2230 
MG 2240 
HG 2260 
HG 226 C 
HG 22 7 C 
HG 2280 
HG 2290 
HG 2 3 0 C 
HG 2 31 -4 
HG 2*20 

HG 2J3C 
HG 2340 
MG 2*50 
MG 236 C 
HG 23 70 
MG 2 3 H C 
HG 2 39 C 
MG 2400 
HG 24 1C 
HG 24 2C 
HU 24 1C 
MG 24 4 C 
HG 2 *♦ SC 
HG 2460 
m G 2 4 7 C 
MG 2 4 8 C 
H G 2 4 9 0 
HG 2 6 0 C 
HG 2 > 1 C 
MG 2G 2 C 
mu ;g3< 

h;i 2 ■ ' 4 t. 
M(» JVu. 

MM 2iM 

MG 2 1 7 < 


on o o o o 
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l 


120 


121 


l?2 


1 23 


C 

C 

C 


W N I IF ( 69 152} *ALLfPHLOSS( J) t ASmMJ) *XHSMM( J) 

UMM 

TT=wALL+P M L055 ( J) ♦XSMM(J) ♦ SUMM 
AA = WALI./TT 
ft H = XSMM ( j ) /T T 
A C - ^ U M m / r T 
Al)sHfiLOSS( J) /TT 
WRITE (6* 151) A A 9 AO « AM t AC 
CONT tNUfc 

PRINT OUT single ION DENSITY PATIO 

WRITE ( 6 * 15 3) NNIO 
SO ( 1 0 • JJ « T I ) =RNI0*XN0 
PRIF=F'R (5*0 
WPITE (6*161) 

PPIO=FR { w« 0) 

PRTHsFR (Oft) 

WRITE ( 6 * 16h) SUM l f PPSUMl fTSU^l *PRIF fTSSfPRlG 

Sf. V = 5f V fSUMlopPSUM l ♦ T SUM 1 *PR I K ♦ TSS 4 >P9 IG 

TT = SUMl ♦ ] SU'4l ♦ TSS 

A A =SU W 1 / T T 

AP-TSU^l/TT 

AC=TSS/TT 

WRITE (6*154) A A » AH t AC 
WRITE (6*162) 

WRITE (6* 167) Wl 05S*SUM?«PRSUM?* TSUM292RIH 

TT = Sijm2*wl 0 SSWSUM2 

A A aWl.OSS/T T 

AP*SUm2/ TT 

AC=TSUM2/TT 

WRITE (S*l46) A A ♦ A m * AC 

CALCULATION OF SlNGlK PL T A /NT UTK At. RATIO 

00 1?3 J=1 *2 
no l?i 1 = 1*2 
no l?l *=1*21 

SX(I»K)=SMl’U,J»ItK) 

SE ( I *K ) =SE l l M ( J, I *K ) 

FAIL T I N T t G (Sr *SF 9 RNIO 9 SUMI f^RSUMl ) 

no 12 ? 1=1*? 
no 122 k = 1 *21 

S x ( I *K ) = SX 1 m? <J9 I 9K) 

SF ( J ,K ) =SF lM? < J« t *K) 

TALI. TINTEO ( 5 x *SE * l • *SUM?*PRSU M ?) 

*(_OSS = VR/OA/F ? 

ra]mo( j) = s u m 1 / ( Sum? ♦ wloss j 

* *< I TE < 69 156) T M J ) 9 R A 1 M 0 ( J ) 
v w I T E ( 6 * 1 6 l ) 

WRlTF (6*163) *UMl,PHSUMl 
SF V = SE V ♦ S 1 J R 1 ORRSUM 1 
WRITE (Of 16?) 

*wP 1 T E (09 lhh) WLOSSfSUM^fPRSUM^ 

T T = WLOSS ♦ SUM? 
ft A - w 1 . 1 ) S S / T T 
Ah^Sijm,V I T 
V.*vl P (6*UM) Aft* AH 
PONT I M)E 

n ( ?f j.j 9 1 n ■= t rnio*ma 1 mo ( n 4 RA i«*’ « ? 1 » *vp*xno*i ou , 

XNm = - ft 1 mO ( 1 ) 4m- ft I MO ( 2 ) 

if ( [ i .(1 r . 1 ) * nm=2. **nm-sp is* .u* 1 1 ) / 100 • 
CALCINATION of » ) 0 1 1 R l f / U 1 1 1 T R A I . RATIO 


HG 

2640 

HG 

2650 

HG 

2660 

HG 

?67Q 

HG 

26AQ 

HG 

2690 

HG 

2700 

HG 

2710 

HG 

2720 

HG 

2730 

HG 

2740 

HG 

2750 

HG 

2760 

HG 

2770 

HG 

2780 

HG 

2790 

HO 

2800 

no- 

2810 

no 

2820 

HG 

2830 

HG 

2840 

HG 

2850 

HG 

2860 

HG 

2870 

HG 

2880 

HG 

2890 

HG 

2900 

HG 

?9 1 0 

HG 

2920 

HG 

2930 

HG 

2940 

HG 

2950 

HG 

2960 

HG 

29 70 

HG 

2380 

HG 

2990 

HG 

3000 

HG 

3010 

HG 

3020 

HG 

3030 

HG 

3040 

HG 

3050 

HG 

3060 

MG 

30 70 

HG 

3080 

HG 

3090 

HG 

3100 

HG 

3110 

HG 

3120 

HG 

3130 

MG 

3140 

HG 

3150 

HG 

1160 

HG 

3170 

HG 

UB0 

Mi 

U 90 

HO 

32 0 0 

HO 

32 10 

HO 

12 2 0 

HG 

12 30 

MO 

12 4 0 

HO 

1 2 5 0 

HO 

12 6 0 

HO 

32 7 0 

H < 1 

3 2 M 0 

M( ' 

1 M0 


CAM Y pit Mi |'iA(V«S l Or’* 1 # • S U M 1 9 P R 1 • O’ ' 1 ) 


n r> n n o n r> 
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r 

it* 


CALL Y t AM ( * (S>. \c . Sr 

CAL! SU^IT (S<1»V»SL 1 » V « /> , k* 1 *0 • T SUMl i8«C> 

WI.OSS - l .4 1 4 ^»Vh/VA/F?» 

RK sfl • 

SD <S* J -it T I ) ** (»A J mp < 1 ) ♦MAI m?j | i» > ) o 1 00 . 

nn 1^4 1=1 *? 

1?4 gr ^ wk^PAO^O ( t ) *SUMM 

HI “0 . 

00 1 ?S 1=1 *NG 
1 2S PL -PL ♦ gp ( I ) 

P?0= ( FljMl ♦ I SUM 1 ♦ HK*Pl ) / tt[ OSS 

PPFniCTinN OF M; XT I lOoPt F ION/NLuTHaL rOP USK TN PPfDlCTION 
OF KPUTPAI. OF NS I T Y 

IF ( I I , FO . 1 ) 00 TO 1^6 

n " T I - 1 

tK '?- (?* »9?fi-sn (^JJ.Il)SODC, 

00 To 1?? 

l?o xm?-^**p?o 

1?7 con T I MJF 

c n ( 1 1 • jjm I ) "P/to*xno 
sr, {? % jj* r I ) --k/M}«*xa 

XI. ( 1 « ..ij * I T ) * VP* 1 • 4 1 A^OXNO^ 100* 

sru o * jj » r i ) = Jr ij* l oo * 

W P I T F nooot K/NFUTWAL w A T I 0 • FTC* 

m p I TP (8*l‘>o> H?0 
(MIT) = P(UM<0 
um-k K (M fC ) 

V P I T K ( h t t m 1 ) 

WPJTF ( 0 • M'S) 


C 

C 

C 

r 

c 

c 


S *■ V “ s *■ .jK 1 ♦ So* ? •” t 

T T “ r 1 l*M ♦ ‘-ij-* ^ ♦ f Si 1 ♦ ■' < ♦ k L 
A A - SUV 1 / T r 
/IfV-SIJM? / I T 
A C = T S U M ! /IT 

a o - p < / t r 

M- iv i / T r 


st ,*a ] ♦ p p s u m l * SUM/* * pr<st»M^ , i sup i « pn • pk » ppsijmm 

4 TSU ■* \ ’* PN ♦ ( u * k ♦ p L ) <} PPSLIMM 


rfF-* I U 
wH | T F 
WO [ T‘ 
WiM I L 
ww rt- 


( G t 1 s 7) 
( 6 • 1 'M ) 
(A» 1 L )M 

( 0 ♦ 1 ^ r 5 ) 

( h • l * h ) 


AA*AM,AC*AI) 
k I t s P S W ►-* f A 

A e 


*1 OSS 

!U- AM CUPPFNT 

xx / x = I *l. < ; ♦ j J* l T) ♦ M. ( 3 * JJ* I I ) <V . ) * 1 .M “ 16 »nPC 
MASS LOSS g.A T e 

To ri ISS-: Ml (l * J J * T I ) <*oos ♦ { «I i ? • JJ* II ) ♦XL ( 3 * J J • I I ) ) r 'OPC) * 

*M - If, 

l; U - UU/t'J| <' V' 

S 1 M ‘ * * J ■ 1 ♦ I P - 1 * M 


1 


f V/ I ON 


M'f (' m. ,o S I 


. I K ;. 


v / ( * * « k / \ , * J - 1 ■' ) ° 1 r,t : v oj_ 


' ■ . ! ‘ i v 

: m 1 * • ■ 

1 *• 1 j 1 


- 1 1 


) A- * I 

1 t 7 


HO 

MQO 

HO 

mo 

HO 

3 

HO 

3 no 

HO 

3 3 a 0 

HO 

3300 

HO 

UGO 

HO 

3 3/0 

HO 

3380 

HO 

3 3 S 0 

HO 

3*00 

HG 

3 1 0 

HG 

1 A ^0 

HG 

3^30 

HO 

.3440 

nO 

3 4 S 0 

HO 

3400 

HG 

34/0 

HO 

3 4 8 0 

HG 

3490 

MO 

3 S 0 0 

HO 

3 S 10 

HO 

30?0 

HG 

JG 30 

HO 

3 ‘j 4 0 

H G 

3 b SO 

HG 

3 SG 0 

MG 

3 *>7 0 

HO 

3 b «0 

HQ 

3^90 

HO 

3 GOO 

HG 

3 G 1 0 

MG 

3 h? 0 

MO 

3 S 30 

HO 

30 4 0 

H G 

3 oS 0 

HG 

3000 

HG 

3 O /0 

Ml, 

30 8 0 

HO 

3 o 9 0 

M G 

3/00 

MO 

. 3/10 

HG 

3 /^U 

MG 

3/30 

HG 

3 / 4 0 

HO 

3 /SO 

MG 

3 /OO 

Mb 

3/70 

MG 

3 /HO 

HO 

i / 9 0 

MO 

3 3 0 0 

Hu 

3 310 

1 Hr, 

H V 0 

MO 

{ < <0 

HO 

G-« 4 0 

HG 

3 .AS (: 

Mii 

< *#xy 

Ml 

- 3 / ) 

O < 

} (•( u 

l‘h 

A ‘ 4 - ‘Mj 

MI, 

i #111, 

H , , 

1 • 1 0 

*M . 

< ) 

Ml, 

< #M , > 
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c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


WP ITE (6*159) XX * x * uTL * FV * ARC I 

CHECK FOP PLASMA NEUTRALITY 

FP*ABS ( <XN M * AN2*RNlfl) /XNl-l • ) 

IF ( IF L At? #E «. • l • A NO, tH #GT. .03) GO TO 108 
PH I NHO = PkI NPfUOP 
IF { IFLAG.NF *0 ) GO TO 138 
12R CONTINUE 

1?9 CONTINUE 

INTERPOLATION POP Tm£ ELECTRON DENSITY WHICH yields oe« 
LOSS RATE 

00 13? I Y = 1 • NS 

LOSS PATF 

no 131 JZ=l*NO 
00 130 J Y - 1 * N S 

PYP ( JY) - ( XL ( 1 * JY t JZ ) *OPN* ( X L t 2 • J Y • J Z ) ♦XL (3 
1 )*l.*F-16 

130 CONTINUE 

WPIIP (6*ln0) (PYP ( JK ) * JK= 1 * NS > 

00 131 JA-lihC 

CALL A I Tk f N (P fP * *NF *NS • 2 • PNO ( JA ) ♦ PNE ( JA * JZ ) ) 
IF (HNF(JAfjZ) • L I • 0 . ) P N F ( J A * J / > = 0 • 

131 CONTINUE 

CALCULATE N* */N + 

00 13? 12—1 • NO 

DPP ( I Y ♦ I 7 ) =SD ( * * IY,I/) / Sl) (AfIYfI7) 

13? CONTINUE 

C A l.L CWSPl. T ( PNC. 4 OPP ♦ A.NE * PNfc * 0 • x 4 \ J * / * f ;!) * NS i SJ C ) 


HG 

3960 

HG 

3970 

HG 

3980 

HG 

3990 

HG 

4000 

HG 

4010 

HG 

4020 

HG 

4030 

HG 

4040 

HG 

4050 

HG 

4060 

HG 

4070 

HG 

4080 

HG 

4090 

HG 

4 100 

HG 

4110 

HG 

4 1 ? 0 

HG 

4130 

HG 

4 140 

HG 

4150 

:hg 

4160 

HG 

4170 

HG 

4 1 8 0 

HG 

4190 

MG 

4 20 0 

HG 

4210 

HG 

4220 

HG 

4230 

HG 

4240 

HG 

4250 

HG 

4260 

HO 

4 ?7 0 

HG 

42«0 

mG 

4 29 0 

HG 

4 3 0 0 

MG 

431 U 


no 13* I * = 1 ♦ 1 3 
on 133 I F - l * R 

133 Y(TF)=TTI(IK*1F) 
no 134 IF =1 *NS 
no 134 IF *1 .h-O 

134 DPP ( IF • IK ) -SO ( l* • IK * IF) 

CAt L C P S P L T (PNO*OPP* XNK.fPNE*G«X .Y.Z*ND*NS*NC) 

135 CONTINUE 

DO 137 L=1*NC 
no 1 3* k = u?* 

PXP(K)=UTT (L»<) /PNO(L) 

13 6 PYP <K) sPDtLf K ) 

TAIL MAP A ( S « P \P *PYH • 1 ♦ ?9 *H|. * HH * VL • V' 4 *9R • NT ♦ / * 1 ) 

137 CALL M A PH ( 5 f pX p . p yp • 1 * ?5 ♦ ML * M H • VI * VM • NR * OT ♦ / ♦ 1 ) 


HG 

HG 

HG 

HG 

HG 

HG 

HG 

HG 

HG 

HG 

HG 

HG 

HG 

HG 

HG 

HG 

HG 

HG 

MG 

HG 

HG 

HG 


138 CONTINUE 

139 FOPmaT (?*• 1GHITKPATION FOR NO 

140 forma r ( 4 p i o * a ) 

1 format (PPio.4) 

14? FORMAT (ISIS) 

143 KORna 1 (AF10.4.F10.P) m ’ 

14 4 FT MM A T (4HV/A* *F*. 3 *9< • .JHtD* «FG. 3* 3HF 1 R *9X • ?MT ^ f F 4 , 1 « 9X « 4 mP /M-= .► ;j ,H(, 

1 ? ) Mfl 
\ 4 G F 0 9 M A I ( / // » IX i 1 <*HVO| 'N4 A i w * 1 1 ^ i i • 4 * 5 X « 1 '* m H ► C uf N (M-3) 3 • ^ \ l , 4 * S H Cj 

1 X t l4MFl.fr Tf HP ( E V ) % f- / • } • 9 X t 16^ ** ' I • f MK *G Y ( E V ) ~ • *■ f • 3 t 5 X ♦ *MNP/MMs * ► MG 

? 6 .*) . H<1 
Us FORMAT (/« l o X • I 1 )h v f;l (C M3) =fE in. 3* 1 '» x ♦ 7HA-NOlT'*F 10. U10X.10MA-* < H G 

1 CM? ) “ • y\ o , 3 • 3 x * 3 hK 1 K ♦ 3 * • 2 = • K 6 • ? ) M(l 


4320 
4330 
4 340 
4350 
4360 
4370 
4 3? 0 
4390 
4400 
44 1 0 
4 4 ? 0 
4430 
4440 
4450 
4 4 6 0 

44 ro 

4 4 8 0 
4490 
4 GO 0 
AGIO 
4 5 2 Q 
4 5 3 0 
4 5 4 0 
4 590 
4*560 

4 s r o 

4580 
4 b 9 0 
40 0 U 
46 10 
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147 KORmaT (///* ?0X* A10*5X*9HNmETA/NO=*F 1Q ,5) 
l4ft FORMAT (1 0* *F 1 0 #5 i 21 X * 2 (F 1 0 .ft * 20X ) > 

1 4 c) FORMAT <//*lOX*23HNhUrRAI. DENSITY (CM-3)* 

10NS*tl4) 

150 FORMAT 
1 *Fft.5> 

15 1 FORMAT 
15? FORMAT 

153 FORMAT 

154 FORMAT 

155 FORMAT 
155 FORMAT 
157 FORMAT 
15? FORMAT 
159 FO»maT 

lev PER 

150 FORMAT 

151 FORmaT 

15? FORMAT ( / ♦ ?0 X ♦ 52HDF NQM I NA TOR TERMS FOLLOWED fiY PRl 
1NY) ) 

153 FORMAT 

154 format 

155 FORMAT 
155 FORMAT 
157 FCRmaT 
15ft FORMAT 


^ t 

( / / * ?0 X * A 1 0 * 5X t rtHNRFS/NO = • F 1 0 *6 * l 0 X * 2 0 HMF AN 
(UX.F10.5*?0Xf3(Fi0.5t?0X)) 

(10X*Fll.4 f 2GX*Fl).4*20X*2(eii.4*lX*E9.2*10; 
(///•? 0 X *5HN*/N0- #F 1 0 . 5) 

( 1 0 X ♦ 3 ( F 10.5*200 ) 

(//♦20X*Ali)*6Xfl UHN*MFT A/NO = * F l 0 • 5 ) 

(///*? 0 ' * 7hn-**/NO=*F 10*5) 

(10X*4(F10.5*21X> ./) 

UOX ,F10.5> 

X * 5 HR E A ♦ F 1 0 » 3 • 3HM A * * 1 C X * 1 2hUT I L I Z A T 1 01 
RF am rON=*F10.5t 10X*5 hIARC=*F10.3) 

{/•2x * jnM/A«y t 3x * p 10*1) ) 

{ / * 20 X * 4QhNUMEkA tor terms (1/SEC) followed 


( 1 0 X * K 1 1 *4* 1 X * E 9 • 2 ) 
(10X*3(F11.4*1X*F9.2*10X) ) 
(10X*4(Fll.4*lX*H9.?*10X)> 

( lOX.F.ll .4*20X*Fll .4*lX*E9.2*lUX> 

(10X*Fll.4*?0X*2(Fll.4*lX*E9.2*10X)) 

<3(F10.4*F10.3)*20X) 


C 

C 

c 

c 


c 

c 

c 


END 

SUBROUTINE CRSPLT (PNO*UPP*XNE*PNE*G*X*H*2*ND*NS*NC) 

THIS SUP POUT I NF INTERPOLATES In thf array DPP TO get 
CONSTANT maSS FLOW RATE (UFT. BY PF’E ) AnD Thr. N PLOTS 


HG 

4520 

HG 

4530 

15HN0. ITERATING 

4540 

HG 

4 650 

REF PATH (CM)srt(i 

4660 

HG 

4670 

HG 

4680 

) HG 

4690 

HG 

4700 

HG 

4710 

HG 

4 720 

HG 

4 730 

HG 

4740 

HG 

4750 

i*F10,6* 10X* 15HHG 

4760 

HG 

4770 

HG 

4780 

r PRI. FRACTIONHG 

4 790 

HG 

4800 

FRACTION (IF AHG 

4810 

HG 

4820 

HG 

4830 

HG 

4640 

HG 

4850 

HG 

4860 

HG 

4870 

HG 

4 8ft 0 

HG 

4890 

HG 

4900 


THE 

THE 


VALUES F 
RESULTS 


COMMON /ft/ UTL(4«?5) ♦Pf)(4f?5) *IK 
DIMENSION DPR(5f5)» PN0(4)* XNF(6)* 

130)* PP ( 30 ) * SR ( 1 0 ♦ 30 ) ♦ X(ft) * Y ( 8 ) * 

interpolation 


PNF (5*10) * 
H ( h ) * / ( 4 ) 


G ( ft ) * PYP ( 3 0 ) • 


C 

C 

C 


DO 102 k ? - 1 *MO 

no 101 kY=1*NS 
101 PYP (KY) sJJPK ( K Y * K 7 ) 

no 102 hx -1 ,mc 

CALL AlT*FN ( XNF * R Y R • MS * 2 * PNE ( KX * K 1 ) * SP ( XX * Ka ) ) 
IE (pmf (KX «K?) *fiT«l 951*1) SP ( KX *KZ ) =0 • 

IF (PNF ( K X * K 2 ) *F Q • 0 . ) SP(KX*K/)=0. 

10? CONTINUE 

FIND max and min 

no 104 jo = i * nc 
no 103 JP=1*ND 

103 PP( JR) -SP( jn* JP) 

104 CALI. MARA ( f- • 5 * PP • 1 * NO * H| * HH * VI, * Vm * X * Y • T I ♦ 1 ) 

PLOT OF DATA 

OAU. MAPA ( 1 *PXP*PYP * l *2H .HL *HH* VL * V! x * Y * T I * 1 ) 
fAll M A PM {1*RXP*PYR*1*25 *kL*HH*VL*Vh«X«Y*TI*1) 

00 112 |A = l* fc) C 

no 105 tMslfNO 

105 PP ( I ft ) =Sw ( I A * I H ) 

CAL! vfNTFRR ( f* • RP *R;i‘ *R YR • * 25 • 3) 

*x -pro r t A ) 


CRS 

10 

CHS 

20 

ICRS 

30 

CRS 

40 

CPS 

50 

CPS 

60 

:ck$ 

70 

CRS 

ftO 

CRS 

90 

CRS 

100 

CRS 

no 

CRS 

1?C 

CRS 

130 

CRS 

140 

CRS 

ISO 

CHS 

160 

CRS 

170 

CRS 

180 

CRS 

190 

CRS 

200 

CRS 

210 

CHS 

220 

CRS 

230 

CRS 

240 

OS 

250 

CHS 

260 

CRS 

270 

OS 

2ft0 

OS 

2 90 

OS 

3u Q 

c«s 

no 

CRS 

120 

os 

3 10 

CRS 

340 

CRS 

35 U 

CRS 

160 
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WPITF (6*113) XX 

CPS 

370 



rfRlTF (6*114) X(l)fHU) 

CRS 

380 



DO 10S lX=lt?S 

CkS 

390 


106 

wRITF (6*11*) PXP ( IX) *PYP ( IX) 

CRS 

400 



WRITF (6*116) (PN'K ( lAf II ) *11 = 1 #N0) 

CRS 

410 



WPI TF (8*117) (G < IX) *PP ( IX) ♦ IX = 1 » NO) 

CRS 

420 

c 



CRS 

430 

c 


PLOT OF UTILIZATION VERSUS DISCHARGf POWER 

CRS 

440 

c 



CRS 

450 



IF ( I**FU.9.0R.I*.FQ.n) GO TO 107 

CRS 

460 



GO TO 111 

CRS 

470 


107 

IF (IK.F0.13) GO TO 109 

CRS 

480 



no 108 L = l *?5 

CRS 

490 


108 

UTL ( 1 1 ti.) =P YP ( L ) 

CR C 

500 



GO TO 111 

c^s 

510 


109 

DO 110 L* 1 *?s 

CRS 

520 


no 

PH ( I A *L ) =PYP (L ) 

CRS 

530 


111 

CONTI NUF 

CRS 

540 



CALL MAPw <?rPXP*PYP*l *?5*HL*HH*VL*VH*A*H*Z.1> 

CRS 

550 


11? 

CALL MAPA ( 2 * PXP • J YP • l • 26 • HL *HH * VL • VH * X * M * Z » 1 ) 

CRS 

560 



CALL MAP A ( 4 *PXP *PYP* l * ?s , rtL *HH* VL * VH* < *H* / * 1 ) 

CRS 

570 



CALL MAPM ('♦ ,PXP*PYP* 1 *?6,hl*HH* VL • Vh* X *H*Z * i ) 

CRS 

580 



PE TURN 

CRS 

590 

c 



CRS 

600 


Ill 

FORMAT (//*?«<» 1 3HP 1 CKF.O V A LUE = * 2 X * E 1 U • 3 * 4mM I LL I A MPS ) 

CRS 

610 


114 

FORMAT ( / * 1 0 X «Al0*20X*AlO) 

CRS 

620 


ns 

FORMAT (]2X*F*.3* 1 6 X * F 1 1.4* 1 0 X t E 1 1 .'♦) 

CRS 

630 


ns 

FORMAT ( / *SX * 1 0 M f LEC.DFNS . • 8 ( ?X * F 10 * 3 > ) 

CRS 

640 


117 

FORMAT (PF10.3) 

CRS 

650 

c 



CRS 

660 



END 

CRS 

670 



SUBROUTINE XINTFRP (X *Y* X I * Y I • N I N * NCUT « InTFRP) 

INT 

10 

c 



INT 

20 

c 


ThIS SiJRWOUTlNF RETURNS -NOUT- POINTS WhICh ARE INTERPOLATED TO 

INT 

30 

c 


ThE -INTEPP- f>F GREF FROM -x,Y- 

INT 

40 

c 



INT 

50 



DImenSTON X ( N I N ) * Y(NIN)* xI(hOUT)* YI(NOUT) 

INT 

60 



nx= ( X (NtN) -x ( I) ) /FLOAT (UnUT-1 ) 

INT 

70 



XX = X ( 1 ) 

INT 

80 



DO 101 1 = 1* NOUT 

INT 

90 



XI ( I ) = XX 

INT 

100 



CALL AITKFN (X*Y*NIN* TNTERP*XX*YY) 

INT 

no 



IF (YY.LT.O.) YY=0. 

INT 

120 



YI ( I ) =Y Y 

INT 

130 


101 

xx-=xx*nx 

INT 

l 4 0 



RETURN 

INT 

150 

c 



INT 

160 



ENO 

INT 

170 



SUP POU TINE YTNTEG ( S I O A * S TGNRG * PkOP * SUM * PR SUM ) 

NTG 

10 

c 



NTG 

20 

c 


evaluation OF REACTION naTF 

NTG 

30 

c 



NTG 

40 



COMMON /A/ NP *NM f T *P^ fNRG 

NTG 

50 



DIMENSION STOMA (2*21 ) * STG^RD(?*?n* SMPlJf SM21) 

NTG 

60 



REAL NP*mm 

NTG 

7 0 

c 



NTG 

80 

c 


PRIMARY F | FO T RONS 

N TG 

NO 

c 



N T G 

100 



no 101 T = 1 * 2 1 

N T 0 

1 1 C 



SX ( n - N !*'*■• A (1*1) 

NTfi 

12C 


101 

SF ( I ) -Slr, f ,v„( 1 , t ) 

N T < i 

1 K 



CALL A I ^ • S X • 2 1 * ^ * PR T NRG * R X ) 

N TG 

1 4 C 



PSUM-rNP^Pt?n*^MX 

NTr, 

1M 

r 



n r g 

1 6i 


oon n n n 


i 


MAX^FLLIAN FLECTIONS 

DO 102 T ~ 1 • 2 1 

S X ( r ) -SIGMA <2.1 ) 

102 SF ( I ) -SIGNRG f 2# I ) 

CALL AITKKN < SF * SX t 2 l f r. . T . TK ) 

TSUM=PPOP*NM*TP 

Sums t Sum ♦ p sum 

PhSUMrPSUM/SUM 

RF TURN 


SUBROUTINE SUM IT ( S I G * SGE • N • PROP * TSUM • SUMS * PPSUMS ) 

TRANSFER routine for excited states 
common /A/ MPtNM, TFMP*PRTNRG 

DIMENSION SIG<N.2*2l>. SGE ( N ♦ 2 . 2 1 ) . PpOP(N)* SUMS(N). 
1 X ( 2 • 2S ) • SE(2*?S> 

REAL NP.NM 

CALCULATES TmF. SUM OF INTERNALS 

TSUM=0. 

DC 102 1=1. N 

no ioi j 2 1 . ? 

DO 101 K = 1 .21 

SX(J.K)aSlRd.J.K) 

101 sf(j.k>=sgf (I.j.k) 

CALL YTNTF6 (SX.SE. PROP(I) iSUMS(I) .WRSUMS(I) ) 

102 TSUMs T SUM ♦SUMS ( I ) 

RE TURN 


SUM 

SUM 

SUM 

SUM 

SUM 

PRSUMS(N). SSUM 
SUM 
SUM 
SUM 
SUM 
SUM 
SUM 
SUM 
SUM 
SUM 
SUM 
SUM 
SUM 
SUM 
SUM 
SUM 
SUM 


SUBROUTINE AITkEN ( X ♦ Y * N ♦ K * XP . YR ) 


A I T k f ' M INTERPOLATION SUBROUTINE 
CALI. TNG SEQUENCE . • . 

CALL AITKFN(X.Y.N.K.XM.Yh) 

X IS A ONE niMFNSTONAl ARRAY Of- l NOEpENOEnT 
VART ARLF ( INCREASING OR DECREASING) 

Y IS A ONE niMFNStONAl ARRAY OK Ob Pr NDENT 
VAPTAGI.F 

N IS MO. OF X.Y PAIRS 

K IS DEGREE OF IMTERPOL AT ING POLYNOMIAL (MAX = 10 1 
<p T c TNOFP. V A R [ ABLE ARGUMENT 
vr IS INTERPOLATED RESULT 

TYPE. O I MF nS I ON AND LABELED COMMON STATEMENTS 

DIMENSION X(N). Y ( N ) • x X ( 1 1 ) . YY(il) 

K 1 2K ♦ 1 

IF ( X ( N ) < ( 1 ) ) 110.101*101 

101 IF < x « - X ( 1 J ) 102.102*103 

102 LL = 0 

00 TO 1 1 Q 

10 1 IF (X(N)-XR) 104.10 A. I 

104 I.L^b-Kl 

or Tr iiR 


AIT 
ATT 
»*A I T 
AIT 
AIT 
AIT 
AIT 
AIT 
A I T 
AIT 
Air 
A IT 
AIT 
AIT 
AIT 
AIT 
° * A I T 
A I f 
A IT 
AIT 
A I T 
AIT 
AIT 

A t r 
A { T 
ATT 
Air 
A IT 
A IT 
A I T 
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RK" 

OW 


105 LL= 1 
LU = N 

106 IF (t.U-LL-1) 117*1 17*107 

107 Ll*(LL*LU)/2 

IF (X(LI)-XS) 108*108*109 

108 LL=L I 

GO TO 10m 

109 LU=tI 

GO TO 106 

110 IF (Xtf-X(l) ) 111*102*102 

111 IF (X(N)-XH) 112*104*104 

11? LL=l 

LU = N 

113 IF (LU-LL-1) 117*117*114 

114 LI*(LL*LU)/2 

IF <X(LI)-xe> 116*116*116 

115 LU=LI 

GO TO 113 

116 LL=L T 

GO TO m 

117 LL -LI - < * l ♦ 1 ) /2 

IF ( L L ) 102*119.118 

118 IF ( LL +K 1 -M ) 119*119*104 

119 DO 120 I “ 1 * K 1 

I 1 S LL ♦ I 

XX ( I) =X ( ID -XR 

120 Y Y ( I ) = Y ( I 1 ) 

DO 1?1 1 = 1** 

DO 121 J = I * K 

IF (XX(j*D.FO.XxtI)) GO TO 122 

121 YY(J*l) = (l./(XX(J*l)-XX(IH) ft (YYU)nX{J*l)-YY(J*l)nX(I)) 
Yh = YY (Kl > 

RF TURN 

12? *R I T F (6*123) (X(t)*I = l*M) 

wPITF (6*123) ( Y (I) *1 = 1* M) 

WRITF (6*124) 

As 1 ./ ( XX ( J* 1 ) -XX ( I ) ) 

CCC=3 • 0 A 
RF TURN 
C 

123 FORMAT (1 0 ( 2 X « F 10.3) ) 

124 FCR^aT ( 10X* 1 7hT90UHLF. U' Aim.N) 

C 

END 


s. coon 

o. 

6.0000 

0. 

7.0000 

0. 

mG-HG ♦ 

8.0000 

0. 

9.0000 

0. 

1 o . o D o n 

0. 

HO-HG ♦ 

1 1 . 0000 

.62AF-U 

o 

C 

O 

•V- 

.164F-13 

1 4 • OuOu 

. 3 3 9 F -1 3MG-rtG ♦ 

1 6.0000 

.ROhF-1 J 

16.0000 

.6 7 3K-1 3 

20 .00*10 

. 83 7F - 1 

23.6000 

. 1 1 1 F - 1 2 

27.0000 

. 1 3 6 F - 1 2 

A 0 • 0 0 0 0 

. 16 IF - 1 2 Hft -**(»♦ 

93.5000 

. 17PF -12 

i 7 . 0 0 0 0 

. 168F-12 

*0,000 0 

. 1 98F - 1 6MG-MG* 

42.(1000 

•205F-12 

4 4.000 (j 

.21U-12 

6 0 . 0 0 A o 

. 630F - 1 2M‘V‘HCi* 

3 • o o 0 0 

.22SP-14 

4.0000 

.6 i if -14 

4 . i o n o 

. 7 92^ - 1 4M(i-Mti ♦ 

4 . 7000 

.10 if -13 

6.0000 

. 12U -1 3 

s. iono 

.IMF-) 

6.7000 

• 1 6 *}F » 1 3 

6.0000 

.hh-n 

6.300 ,) 

• 2 1 4 F - 1 OiO-hO ♦ 

6.7000 

.244F-1 1 

7 . 0 0 0 0 

, if 6 /U - l 1 

7 . 6 o n i 

. U)H> - 1 

iHG-HG * 

8.0000 

. 3*Hf -1 3 

« . S 0 0 0 

. * M O r - 1 1 

. fMJ 0 U 

, -.2^ - 1 

iM(i-Mf) ♦ 

9.^000 

.460t - ) j 

1 U . 0 0 0 0 

. 6 U 9 r - 1 i 

1 0 ,‘i Jf'O 

, 6 4 F - 1 

ihij-Mli ♦ 

11.0000 

.6891- - 1 3 

12.0000 

. 6 6 * f - 1 A 

1 i . o o n o 

. 7 4 (j F - l 

ihlj-MG* 

6 . 0 o 0 0 

0. 

6 . 00'H) 

0. 

7 . o oo o 

0 . 

Nfj-M(* ♦ ♦ 

8,0000 

0. 

9 . 0 0 0 0 

0. 

1 o . n o 0 0 

f ). 

Mi* - Hi* ♦ ♦ 


AIT 

310 

AIT 

320 

AIT 

330 

AIT 

340 

AIT 

350 

AIT 

360 

AIT 

370 

AIT 

380 

AIT 

390 

AIT 

400 

AIT 

410 

AIT 

420 

AIT 

430 

AIT 

440 

AIT 

450 

AIT 

460 

AIT 

470 - 

AIT 

480 

AIT 

490 

AIT 

500 

AIT 

510 

AIT 

520 

AIT 

530 

AIT 

5 A 0 

AIT 

550 

AIT 

560 

AIT 

570 

AIT 

580 

AIT 

590 

AIT 

600 

AIT 

610 

AIT 

620 

AIT 

630 

AIT 

640 

AIT 

650 

AIT 

660 

AIT 

670 

AIT 

680 

AIT 

690 

AIT 

700 

AIT 

710 

AIT 

720 

AIT 

7 30 

AIT 

740 

AIT 

750 

AIT 

760 

AIT 

770 

AIT 

780 

AIT 

790 


1 

2 

3 

4 
6 
6 
7 
1 
2 

3 

4 
*> 
6 
7 
1 

2 
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5. 

10 . 

16 . 

20 . 

25. 

33. 

42. 


11.0000 0. 

12.0000 

0 . 

14.0000 

16.0000 0. 

1 H. 00 00 

0 . 

20.0000 

23.5000 0. 

27.0000 

0 . 

30.0000 

33.5000 .100E-U 

37.0000 

. J06F-14 

40.0000 

42.0000 ,536f-l4 

44.0000 

.64 8F - 1 4 

50,0000 

3.0000 . 322F- 1 « 

4.0000 

.*♦29^-17 

4.3000 

4.7000 .139E-16 

5.0000 

. 2 1 OF- 1 6 

5.3000 

5.7000 .470E-16 

6.0000 

• 629F -1 6 

6.3000 

6.7000 .1UE-15 

7.0000 

.142F-1G 

7.5000 

fi.0000 .26Sf-is 

8.5000 

.346F-15 

9 . 0 0 0 0 

9. 5000 , C '46F -IS 

10.0000 

• 665F - 1 8 

10.5000 

11.0000 . 942E-1 S 

12.0000 

• 1 27F-1 4 

13.0000 

-3. E-10 

7. 

-2.5 F-30 8. 

. 6 


0. 

0. 

.617 E 
.h34F 
.103F 

• 74 3F 

• 30 4F 

• 522F 

• 197F 

.43 9F 
. 797F 

. 164F 


HG-HG** 
HG-HG* ♦ 
-15HG-HG** 
- 1 4HG-HG* * 
1 3HG-HG* ♦ 
-17MG-MG*>4 
• 1 6MG-HG ♦ ♦ 
16HG-HG** 
-15MG-HG** 
1 5HG-HG* ♦ 
• 1 5H6-HG ♦ ♦ 
-1 4HG-HG** 


-1.5 F-°Q 12. 
-.5 F-30 15.6 

.443 F-13 20.5 
.657 F-13 ?7*5 

.94] F - 1 3 37, 
1.15 F-l 3 


• 1 • 25E-30 14. 

0. F-30 19.25 

.47 E-l 3 22. 
.766 F-l 3 30. 
1.02 F-13 40. 


— 2 . 
- 1 . 

.162 

.529 

.654 

1.1 


3-.0000 

.267F-15 

4,0000 

* • r. ■* r. - i .■) 

• 12 7t-l 4 

. 1 

4.3000 

4.7000 

• ? c; 7F-l4 

5 . 0 u 0 0 

.327E-14 

5.3000 

5*7000 

•521 F-l 4 

6 . 0 0 0 u 

• 6 1 6 E - 1 4 

6.3000 

6.7000 

.P60F-14 

7.0000 

• 9 7 3 F - 1 4 

7.8000 

8.0000 

.13HF-1 3 

8. 5000 

• 1 5 9 1 - 1 3 

9.0000 

9.8000 

• 20 3F -1 3 

10 .O 000 

.226F-1 3 

10.5000 

11.0000 

• 2 7 2 F - 1 j 

12.0000 

•318F-13 

13.0000 

5.0000 

.988F-14 

6.0000 

•121F-13 

7. 0000 

8.0000 

• 9 0 7F- 1 4 

9.0000 

. 7 44F - 1 4 

10.0000 

n .oooo 

* 4 4 1 F - 1 4 

12.0000 

. 3 33E-U 

14.0000 

16.0000 

•191F-14 

18.0000 

• 123F-14 

2 0 • 0 u 0 J 

23.5000 

• F 27f - 1 5 

27.0000 

. 72UF-15 

30 . 0 0 0 0 

33.5000 

• 4 2 3 F - 1 5 

37.0000 

.286F-15 

40.0000 

42.0000 

0 . 

44.0000 

0 . 

5 0.0i)00 0 

3.0000 

.31 9F-14 

4,0000 

. 3 76F - 1 4 

4.3000 

4.7000 

. 3 9 6 F - 1 4 

5.0000 

. 3 9 8 F - 1 4 

5 . j 0 0 0 

5.7000 

•393F-14 

6.0000 

.396F-U 

6 . 3 0 0 u 

6.7000 

. 388F - 1 4 

7,0000 

. 3 6 4 F - 1 4 

7.5000 

8.0000 

• 367F - 1 4 

6.5000 

. iSHF-14 

9.0000 

9.5000 

•339F-14 

10.0000 

• 3 10E-1 4 

10.5000 

1 1 .0000 

• 3 1 2F - 1 4 

1 2 . 0 0 0 u 

,^94F-l4 

13.0000 

5.0000 

0. 

6.0000 

• 4 3 Y fe - 1 3 

7.0000 

8.0000 

• 4 04F - 1 3 

9.0000 

. i29F-n 

10.0000 

1 1 .0000 

.21 2F -13 

1 2 . 0 0 0 0 

. 165F-13 

U.OOOO 

16.0000 

. 106F-13 

16.0000 

. 7 9 9 F - 1 4 

20.0000 

23.5000 

.433F -14 

27.0000 

.312F-14 

30.0000 

33.5000 

•168F-14 

37.0000 

. 1 44fc -14 

40.0000 

42.0000 

.1 13F-14 

44,0000 

.P16F-15 

50.0000 

3.0000 

.96 IF- 1 4 

4 . 0 0 0 u 

.125F-13 

4 . 3 o 0 0 

4 .7000 

• 1 36F- 1 .1 

8.0000 

. 1 3 8 F — 1 3 

5.300i) 

5*7000 

. 142F-11 

6. 0000 

•1432-13 

5.3000 

6.7000 

.143F-13 

7.0000 

» 1 4 2 K • 1 3 

7.500 0 

8*0000 

. 135F-1 3 

8,4000 

. 137F-13 

9.0000 

9.5000 

• 1 32F - 1 3 

1 0 . 0 0 o u 

. 129F-13 

10.5000 

1 1 .0000 

• 1 2 IF - 1 3 

12.0000 

. 1 1H* -1 j 

1 3 . 0 •) 0 0 

5 .0000 

.948F-14 

6 . 0 0 0 0 

. ^6 IF -14 

7.0000 

8.0000 

• 4 1 0 F - l 4 

9 . 0 0 C 0 

• 6 ■ J 4 r - 1 ** 

1 0 • 0 J o o 

11.0000 

. 7 0 8 F - 1 4 

12.00 0 0 

.606F-U 

1 4 . n no o 

1 * , uono 

• 6 4 0 F - 1 4 

1 " . 0 U 0 0 

. 6 '♦ 1 F - 1 

2 '>.0 0 »o 

23.5000 

.^96F - 14 

27.0000 

*527 F-l 4 * 

in.oooo 

33.6000 

• 3 1 9F - 1 h 

3 7.0 J00 

. 1 6 1 F - 1 4 

4O.CQO0 0 

42.0000 

0 . 

44.0000 

0. 

50.0000 0 

1.0000 

• 2 1 4F - 1 4 

<* . 0 00 () 

• 3 1 8 ► - 1 4 

4 . 1 i ) 0 0 

4.7000 

. 1 c 6f -14 

4 • r. u o o 

. 3 7 i t - 1 4 

5.3000 

6 . 70. JO 

. 198F-14 

4 . OoOi* 

. 4 0HK-1 4 

6 . 3 iJ '• J 

6,7000 

.426F-1- 

7 . 0 00 0 

• .Vr -U 

7.^000 

» . Poon 

• 4 4 6 F - 1 <* 

* , 5 0 0 0 

• *♦ 4 '* h - 1 #« 

4 , «i , » r. i j 


F 
F 
F 

F 

F> 

F 

1.33 F 
. 177F 

• 4 0 5F 
. 7 1 7F 
.1 17F 
. 1H1F 
, 24 9F 
.364F 

• 1 0 7F 

• 5 76F 
.22 HF 
. 103F 
. 5 7.9 F 
» 1 H8F 

.3672 

• 399F 
» 39 3F 

• 3 7 6 F 
.34MF 

• 32 1 F 

• 2 7 8 F 
.482F- 

• 2 h 0 F - 
. 136F- 
. 655F - 
.22 7F- 
. 1 3 1 F 

. J65F 
. 130F- 
.UOF- 

• 1 4 3F 

• 1 4 1 F - 
.134K- 

• 1 26 F - 
. 1 1 2r - 

• 4 h 2 F - 
.7 1 9 F - 
.626** - 
.6 43F- 
, 4 4 4 F - 


- 3 OHG ♦ 
3 0 H G ♦ 
l 3HG ♦ 
1 3hG* 
- 1 3 m G * 
- 1 3HG ♦ 
1 3hG* 
1 4HG ♦ 
1 4hG* 
14MG* 
1 3HG ♦ 
1 3HG* 
1 3rlG* 
1 3hG* 


HG* ♦ 
HG ♦ ♦ 
HG** 
HG** 
HG* ♦ 
HG ♦ ♦ 
HG ♦ ♦ 
HG ♦ ♦ 
HG ♦ ♦ 
HG* ♦ 
HG** 
HG** 
HG** 
HG* ♦ 


.JJ4F- 

• 3 H 4 F - 

• 4 1 6 F - 

4* - 

• 4 S 1 f — 


•13HG-HGM 3PQ 

• 1 4HG-HGM 3P0 

• 1 4 H G - H G M 3P0 

• 1 4HG-HGM 3P0 
■15HG-HGM 3p0 
1 1 8HG-HGM 3P0 

HG-HGM 3P0 
14HG-HGM 3P0 
1 4HG-HGM 3P0 
14mG-HGM 3P0 
14HG-HGM 3P0 
14H6-6GM 3P0 
14MG-MGM 3P0 
) 4HG-MGM 3P0 
13HG-HGM 3P? 
lJHG-HGM 3P2 
1 iHG-HGM 3P2 
1 4 H G - H C M 3p<? 

l4H(3-h6M 3P2 
1 4 H(>-HG m 3P2 
3P2 

13HG-H6M 3P2 
1 3HG-HGM 3P2 
13HG-HGM 3P2 
1 3 m G-mGM JP2 
1 5HG-HGM 3P2 
1 3hG-mGM 3P2 
13HG-HGP 3P2 

1 **h(i-MPW 

1 4H(i-mGW 
1 4HG-HGP 
i JP 1 

l'ihG-HGW 3P J 
HfWlOP 3P ) 
3P 1 

14 m^-HGP 3P | 

1 •♦MG-,109 3P 1 
1 

1 *♦ Mh- Mu 9 
1 ♦ MG - . 1(» V 


3P1 

3Pl 

3P1 


3PJ 

IP) 

*Pl 


1 

2 

3 

4 

5 

6 
7 
1 
2 

3 

4 

5 

6 
7 
1 
2 

3 

4 

5 

6 
7 
1 
2 

3 

4 

5 

6 
7 
1 
2 

3 

4 

5 

6 
7 
1 
2 

3 

4 

5 

6 
1 

1 

z 

3 

4 
6 






r~rrrr. 
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9. soon 

.452F-14 

10.0000 

.461F-14 

1 0 • b 0 0 0 

11 • 0 0 0 0 

• 4 4 8 F - 1 4 

12.0000 

» 4 4 1 F ™ 1 4 

13.0000 

s.0000 

0. 

6 . U J ( J U 

o • 

7.0000 

M . 0000 

• 3 4 9 F - 1 3 

9,0000 

, 6 7 9 F - 1 3 

10.0000 

1 1 . 0000 

. 106F-1? 

12.0000 

. 130F-12 

14.0000 

16.0000 

•179F-12 

1 H , 0 0 0 0 

, IH9F-12 

20.0000 

23. SOOn 

•209K-12 

27.0000 

.216F-12 

30.0000 

33. bOOO 

• 2 2 3 F - 1 2 

37.0000 

.224F-12 

4 0 , o 0 0 0 

42.0000 

.224F-12 

44,0000 

.224F-12 

b 0 * 0 0 0 0 

3.0000 

. 166F-1 3 

4,0000 

.323F-13 

4.3000 

4,7000 

• 4 3 7 F - 1 3 

6,0000 

.485F-13 

5.3000 

5.7000 

• 5 9 4 F • 1 3 

6 . 0 U 0 0 

.9 38F-13 

6.3000 

6.7000 

•737F-1 3 

7.0000 

.777F-13 

7.5000 

p.ooo: 

• 9 0 1 F - 1 3 

8.5000 

.9b7F-i3 

9.0000 

9.5000 

• 1 0 6F - 1 2 

10.0000 

. illF-12 

lU.bOOO 

11.0000 

•1 19F-12 

12.0000 

.127F-12 

13.0000 

5.0000 

0. 

6.0000 

.96H-15 

7.0000 

8 . 0 C 0 0 

• 2 3 3 F - l 3 

9.0 000 

. ibbF - 1 3 

10.0000 

11 .0000 

.590P -13 

12.0000 

. 7 0 1 1 - 1 3 

14.0000 

1 6 . o 0 0 0 

.1 10F-12 

1 8 . 0 0 O 0 

. 1 3 0 F - 1 ? 

20,0000 

23.5000 

.205F-12 

27.0000 

.247L-12 

3 C . 0 0 0 0 

33.5000 

.309F-12 

37.0000 

. 335F-12 

40.0000 

42.0000 

.365F-12 

44.0000 

• 3 7 6 K - 1 2 

SO .0000 

3.0000 

. 1 0 3 F - 1 3 

4.0000 

. 2 1 Oh “ 1 3 

4 . 3 U 0 0 

4.7000 

.29HF-13 

5.0000 

.3388-13 

b • 3 0 0 0 

5.7000 

.433F-1 J 

6.0000 

.4 7 6 F - 1 3 

6.3000 

6.7000 

• 5 7 5 F - 1 3 

7.0000 

* 6 l *3 F - 1 3 

7.60 0 tj 

8.0000 

• 7*2F - 1 J 

8 . 6 0 0 0 

.8 5 4 F “ 1 3 

9.0000 

9 . 5 0 0 0 

.974F-13 

1 o. ooco 

•104F-12 

10.6000 

1 1 . 0 0 0 O 

. 1 1 8F-12 

1 2 . 0 0 0 0 

. 1 UF-12 

13.0 0 0 0 

5,0000 

*682 F -16 

6.0000 

.96HF-1 4 

7.0000 

8.0000 

•484F-1 3 

9. 0000 

♦ 9 6 2 F - 1 3 

1 0 . 0 0 0 0 

it .ooon 

, 0 2 9 F - 1 3 

1 2 . 0 U 0 0 

.1067-12 

14.0000 

1*.0G00 

. 1 49F -12 

18.0000 

. 169K-12 

20 ,0000 

23.5000 

,244t -12 

27 , 0000 

.2868-12 

30.0000 

33.5000 

.346F -12 

3 7 . 0 0 0 0 

. 17 l F - 1 2 

4O.0000 

42.0000 

. 4 0 0 F - 1 2 

44 , 0 jO 0 

.4106-12 

6<) • 0 0 u 0 

3.0000 

. 17mF -1 3 

4.0000 

. 3 2 7 K - 1 3 

4,3000 

4 . 7 0 0 0 

.442K-13 

6.0000 

.4926-13 

6 . 3 U 0 0 

5.7000 

.F 1 OF - 1 3 

6.0000 

.6611-13 

6 . 5 on o 

5,7000 

.7 7>n -13 

7 .000 0 

.6287 - 1 i 

7 • 6 1) 0 u 

8 , 0 0 0 0 

.9Q2F -1 J 

h . S 0 0 0 

.1077-12 

9.00 0 0 

9.5000 

• 1 2 5 F - 1 2 

1 n.OOO') 

. I .10h-12 

10.6000 

11.0000 

.145F-12 

1 2.0000 

• l-v»7 -12 

] 3.0000 

5.0000 

0. 

6.0000 

0. 

7.0000 

8.0000 

0. 

9 , 0 0 n 0 

0 . 

10.0000 

11.0000 

0. 

12.0000 

0. 

1 4 , o ‘JO J 

16.0000 

0. 

18.0000 

0. 

2 0 .00 0 0 

23.5000 

. 722K-14 

2 7 . no no 

.2.Ihh -13 

3 ° .0000 

33.5000 

• 6 7 6 F - 1 1 

3 7 .000 0 

. 7 4 4 7 - 1 J 

a C. . 0 -l 0 0 

4 r ,0 000 


44.0000 

. 1047-12 

^ r . o rl r, (J 

3.0000 

. 570h -16 

4 , o i) 0 0 

.2617-16 

4 , it! f; J 

4.7000 

, ^ a 26 - 1 -> 

6.0000 

. 5 797 -lb 

s . < f ) t) 0 

8.7000 

. 1 6 2 f - 1 4 

6 . 0 0 0 0 

,802V -1 4 

, 9 ) J» 0 

6.7000 

. U5K-U 

7 . IJU8 1I 

• 1 /r't - 1 4 

7 .50 li ) 

8 . o o o n 

• - 1 

8 ,1hPli 

• 7 8 2 F - 1 4 

.O'! " 0 

9.5000 

.ion-} i 

1 0 . 0 00 0 

.11*7-1 H 

i n , ‘ • j 5 ■) 

11.0000 

. 1 49f - l i 

12,0000 

• l -1 3 

1 ',(<■} ti i 

5 . 0 O 0 0 

0. 

r , u ikh) 

0. 

7 . 0 .)0 0 

8.0000 

0. 

9.00 0 o 

0 . 

1 " .6 

ll.oooo 

0. 

1 8 . oo no 

0. 

} u , 0 J -' (i 

l 6 . u 0 o n 

0. 

l M , 0 'too 

(i . 

8i. . fi i) : J 

21.500 o 

• 7 2 2 K - 1 *♦ 

2 7.00 0 0 

, ' <07-1 1 

. r ■■‘O 

7 5.500 f) 

. r 76F - 1 5 

H 7 . () 0 0 .1 

. - 1 1 

<! 

4 2.0 O O 0 

/G'O -1 i 

‘4 «« , 0 0 ■-) t , 

• 1 ' «* * -18 

« , ' » i i) 

7 . 0 0 0 0 

. 7 7 OF - l 6 

•* . •') »'■' ■> 

’’ 1 ; - 1 < 



• 4 6 0 F * 
.4 ;i3F- 
.1 32F- 
, H 1 4 F - 
. 162F- 
. 199F- 
.220F- 

, 2 24 F • 
.222F- 
, 3 72F * 
. b 3 2 F - 
.6H2F- 
.H41F- 

• 1 0 1 F - 
•llbK- 

• 1 3 3 F - 

• 8 6 b F - 
.477F- 
•916F- 

• 1 5 7 K - 

• 2 7 h 17 - 

• Jb 4 P- 
.400 c - 
.247F- 
. J78F- 

• 5 1 rt F - 
.690F- 
.904F- 

• 1UF- 

• 1 4 37 - 

• 3 1 2 F - 

. 796F- 

. 12VF- 

.19/7- 
. JlfaF- 
. 390F- 
. 4 3 3 F - 
. 3 7f>K- 

• 5427 - 

• 7 1 1 7 ' - 
. 9 n^ - 
.1157- 
. 1 38 * - 
. 1 727 - 


3P1 
3P1 
1P1 
1P1 
1P1 
lPl 
1PI 
1P1 
1P1 
IPX 
1P1 
1P1 
1P1 
1P1 
IP1 
1 PI 


14HG-HGP 
1 4HG-HGR 
1 3PG-HGR 
1 3hG-HGR 
1 

•12HG-HGW 
IPhG-mGW 
12HG-HGW 
12mG-HGR 
1 3hG-hGP 
1 3MG-HGR 
1 3HM-MGW 
1 3HG-HGR 
1 2HG-HGR 
12 HG-HGR 
l 2MG-HG R 
l 4HGM-HG* 0 
1 JnGM-HG* 3P0 
1 3HGM-HG^ 3P0 
12 mGM-HG 4 3 P 0 
1 2HGM-HG ♦ 3P0 
1 2mGm-hG ♦ 3P0 
1 PhGm-hG^ 3P0 
1 3HGM-H6 ♦ 3°0 
1 3HG M -H(i ♦ 3P0 
1 3MGM-Hii*3PQ 
] 3 P (i M - H G ♦ 3 p 0 
] iMliM -hG ♦ 3PU 
12HGM.-HG*3P0 
1 2MGM-HG 4 3PQ 
1 3 HGM-HG^ 3 Pr > 
1 3MGP-HG 4 3P2 
1 2HGM-HG ♦ 3P8 
1 2MGM-MG * 3PP 
1 2MoP-MG * 3P? 
1 ft HGM-HG ♦ 3P2 
1 2HGM-HG ♦ 3P2 
1 JmGM-hG* 3P? 
1 3HGM-HG ♦ 1P2 
1 JHGM-MG* 3P2 
1 3Pg 

1 2MGM-MGOP2 
1 2MGM-HG* 1P2 
1 2mgP-ho* 3P2 
HGM ♦ HG 4 *b 


.6267 
. r*b‘- 
.6 7 7 * 
. 1 24F 
. 59 77 
• 1 U>7 

4 2 4 7 V 

.*♦77- 
. -U Itr 
. 1 1 1 »- 
.2217 


. 8‘*8 
. -i 7 7 8 

• 1 8 4 > 

. V * ft 


H GM ♦ 

H(')M ♦ 

■ 1 ‘ihiiM ♦ 

• l 3 m*M ♦ 

• t G U t M ♦ 

• 1 2 HGM ♦ 

■ ] b no M ♦ 

- 1 H M G ♦ 

• | 4 MGM ♦ 

• 1 '♦ M G M ♦ 

• 1 4 mGm . 

- 1 UM * 14 • 

“ 1 <'«(■"♦ M t « ♦ ♦ 1 J 

MG>* ♦ -MG ♦ ♦ 3 

Mi t H ♦ -Hf< 4 ♦ 3 
I'lMli'U-MhM J 
J j i « ( , M » - M i , » • 3 
] \ m t , ** 4 - m-< ♦ ♦ i 

i , • : if, • 4 -mg 4 4 < 

i » S 


MG ♦ ♦b 
MC* ♦ ♦ 5 

MG ♦ ♦ b 
MG ♦ ♦ b 
MG ♦ ♦ b 
MG ♦ ♦ > 
mg 4 ♦ b 

M(, 4 4 b 
MG * ♦ b 

M(l 4 4 b 

MG ♦ ♦S 
MU ♦ ♦ b 


6 

7 

1 

2 

3 

4 

5 

6 
7 
1 
2 

3 

4 

5 
(■ 
7 
1 
2 

3 

4 

b 

6 
7 
1 
2 

3 

4 
b 
h 
7 
1 
2 

3 

4 

b 

h 

7 

1 

2 

3 

4 


in < 
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4.70 00 

5.70 00 
ft, 7000 
A. 0000 

9.5000 

1 1 . 0000 

5.0000 

8.0000 
n.oooo 

lft, 0000 

23.5000 

33.5000 
4?, 0000 

3.0000 

4.7000 

5.7000 

6.7000 
8 . 0 0 0 0 

9.5000 

11,0000 

5.0000 

A. 0000 

1 1 . 0000 
lft . 0000 
2 3 . 5 0 0 0 

13.5000 
4?, 0000 

3.0000 

4.7000 

5.7000 
ft. 7 000 

5.0000 

9.5000 

11.0000 

5.0000 

8.0000 
1 1 . 0000 
lft . 0 000 

23.5000 

13.5000 
4 2 , 0 0 0 0 

3.0000 

4.7000 

5. 7000 
ft . 7000 

5.0000 

9.5000 

11.0000 

5.0000 

5.0000 
1 1 . 0000 
l ft . ooon 
2 3 . 6 0 0 0 
13.50U0 

42.0000 

3 .0000 

4.7000 

5.7000 
ft . 7 0 0 0 

5.0000 

9.5000 
1 1 . >1 0 0 0 

1 1 

1 . 

.0 14 


.6426 -] 5 
. 1 6 2 6' - 1 4 

• 3 1 5F - 1 4 
.5046"- 14 

• 1 0 1 F - 1 3 
. 149F-13 

0. 

• 2 8 9 F - 1 3 
.6706-13 
.119F-12 
.2146 -1? 
.31«F-12 
.374P-1H? 

• 1 1 9 K - 1 3 
.331^-13 
.475F-1 3 

• 623F.-13 
.5175-13 
. 103F-12 
. 1?46 -12 

0. 

, ft ?5f -14 
•339F-13 
.793^-13 
. 1 7 3 r - 1 2 
•279F-12 
. 3 3ft F - 1 ? 
.55SF- 1 4 
. ? 0 1 F - i 3 

• 31 IF-] 3 
. 4 3 1 F - 13 
.595F-13 
» 7 8 ft * - 1 3 
.974 F -13 
. 174F-1 3 
.2) Oh'-Ui 

• 2 3 5 F - 1 3 
,2166-13 
. 1 37F-13 
.51 7 f - 1 4 
,^4F-N 
.8066-14 
. 1226-13 
. 1 3 5 F - 1 3 
. 1 43 c -l 3 
. 1 4 7K-1 3 
.1486-1 3 

• 1 4 ft F - 1 3 
0. 

• lft4F-l 1 
.] Am ft - 1 3 
. 1 Hfif - 1 3 
. 1 1 26-1 3 
, ft ft 4 F - 1 4 

• *♦ 1 - 1 4 

.41 1 6 - 1 4 
. 75..F - 14 

• M 8 26 - 1 4 
a 9 ft fV — 1 4 

. 1 o V -1 3 
. iOftf- -1 3 
.10??-! 3 
ft 1 

1 . 

4.5 


5.0000 

ft . 0000 

7 . o o n «j 

8.5000 
10.0000 
1 2 • 0 0 0 0 

ft. 0 0 00 

9 . 0 0 0 0 
12.0000 
lM.OUOO 

27.0000 

3 7.0 0 0 0 

44.0000 

4.0000 
b . 0 0 0 0 

6.0000 

7.0000 
8 . 5 0 o o 

1 0 . 0000 

12.0000 
b . 0 0 0 0 

9.00 0 0 
1 2 . 0 0 0 0 
18.0000 

27.00 00 
j 7 , 0 0 0 0 

44.000 0 

4.0000 

5.0000 

6.00 0 0 

7.0000 

8.5000 

10.000 0 

12.0000 
6 • 0 0 0 0 

9 . 0 U 0 U 
1 2 , 0 0 0 ti 

1 M . 0 i) 0 0 

2 7 , 0 u 0 0 

3 7.0000 
4 4 . v ) 0 0 0 

4 , no n o 

5 . 0')0 0 
ft . 0 0 0 0 

7 . 0000 

5 .5000 
1 o . 0 Ij f) 0 

12.0000 

6.0000 

9 . 0 0 0 0 
1 2 • 0 0 0 0 
18.0 )*>0 
27.no,M, 
17,0 o 0 <j 
4 4 ,000 0 

4.00 n i, 

u • mu '. >; 

ft . :» 0 0 0 
7 . «) .}0u 

5 . *>(*0 M 

1 (’ . -MHO 
1 • U n *’ j 



.57 OF -15 
.2026-14 
. 3 7 2 F - 1 4 
. 7 2 2 1* - 1 4 
. 1 1 ftF-1 3 
. 1 Mbh - 1 3 
.22 IF -14 

• 4 2 3 F - l 3 
.78/6-1 i 
. 1 '♦ 0 *" -12 

• 2 5 ft K - l 2 
» 3 4 4 F - l 2 
.3846 -12 
. 2 3 7 K - 1 3 
. 3 7 3ft — 1 i 
.51 OF- l 3 
.6586 -1 3 
, o 9 0 K - 1 3 
. 1 11 1: - 1 2 
.1376-12 

0. 

. 1 6 1 F - 1 3 
.4356-13 
.93^6-13 

• 2 1 5 F - 1 2 

• 3 0 5 F - 1 2 
. 34 76-12 

• 1 3 4 K - 1 3 
. 2 3 3 K - 1 3 

• 3 4 ft K - 1 3 
.4686-13 
.6596 -13 
. 8 496 - 1 8 

• 1 0 96 - 1 cr 
.1386 -lj 
.2 196.-1 3 
.2411-1 3 

, 1901.-1 J 
« 1 1 3 K - 1 1 
.6 v / F - 1 *♦ 

.521 F - 1 4 
. lOo) -1 3 
.1866-1 3 
. 1 386-1 3 
. 1 4 4 F - 1 3 

• l 4 h t - 1 3 
. 1 4 h f - 1 \ 

, 1 4 4 F ** 1 3 

[). 

.1 736-1 3 

. 1 9 ft f -1 3 
. X ^ 7t -1 3 
. 988 k - 1 

• 6 6'n - ! . 

. » 8 1 | - 1 <* 


• < g - ; 


.7 . 7 F — 1 

4 

.9|1|-| 


• 9 > ■ 1 ' f - ] 


. 1 9 * 1 

. 1 0 -5 . - ] 

1 

J 

.1 ' ■ — 1 

1 


5.3000 
ft. 3000 

7.5000 

9.0000 
10,500 0 

13.0000 

7.0000 
1 o . n o n o 

1 4 , 0 lj n (, 

20.0000 

3 0 • 0 0 0 U 

4 0,0 0 0 Q 

5 0.0 0 o 0 

4.3000 

5. 3000 
ft. 3000 

7.5000 

9 . 00 00 

10.50 0 0 
1 3.000 •) 

7.0000 

10.0000 
1 4 , (j o o 0 

20 .0000 

30.000 0 

4 0,0 0 0 0 

5 0,0 0 0 0 

4 . 3000 
5. 30 OQ 

6.30 n o 

7.90 0 0 

9 . 0 0 0 J 

10.500 0 

13.0 0 0 0 

7 . 0 0 0 0 

1 0 . 0 0 0 0 

14.0000 
2 0,0 0 0 0 

3 n . o o o ) 

4 0 , 0 00 J 

5 0 , noon 
4 . 3 ■ j 0 0 
- , 3 .) 0 0 
ft ,5000 
7 . 5 P n u 

9 .000 0 
1 0.50 0 0 
1 I.POOu 

7 . n o o o 

10.0 f) n o 
1 4 , <» 0 O I) 
2t I , 0 .'Ml ti 

3 0 , 0 0 »u 
<♦ 0 . fi ■')( ■» 
50.U 0 0 0 
4 . 1 10 0 
. 4 )" 0 
ft , \ 0 ) .! 

7 . ‘ .1 0 •' 

'* . > 5 Oil 4 

1 ■: • u : n m 

• 9 t>:n ) 

i. 1 

* . '« t J 


• 1 1 ftF- 1 4HGM4-HG4 *3 2 

• 24 7F - 1 4hftM ♦ -Hi) ♦ ♦ 3 3 

• 4 / /F- 1 4HfiM*-H6* #.J 4 

,HfeUF->l4H(iM*-HG**3 5 
« 13 38-1 iMftM ♦ -HG^ ♦ 3 6 

• 22 1 F - 1 jHftM ♦ -MG ♦ ♦ 3 7 

. 129K-1 3 h ft 6< - h G ♦ 3 P 1 1 

. 6538 -1 Jh(>W-hG#3Pl 2 
. 1 0 1 6' - 1 2 M ft P - H G ♦ 3 P 1 3 

• 16 7F -12hGM-rtG».3Pl 4 

.285 6 - 1 2hftH-riG ♦ 3P 1 5 

. 3ft 16" - 1 2MGP-MG + 3P 1 6 

.4096 - 1 2MGW-HG ♦ 3P 1 7 

.2/6F-1 3h(iP-HG*3Pt 1 
. 4 1 6 f — 1 3 N ft P - H G ♦ 3 P 1 2 

. 863F -1 3HGW-HG.3P1 3 
. 74 3F-1 3MGP-HG.3P1 4 

. 96 3ft- 1 3HftP-HCW3Pl 5 

.11 76-1 2h(-P -MG ♦ 3 P 1 6 

• 1 50F -1 2HG8-MG* 3P1 7 

. 793h -1 5H09-HG* 1P1 1 

.246F-1 3HGR-HG* IP! £ 

• 6256 -1 3 H ft W - 1 1 fi ♦ 1 P 1 3 

• 1 206 - 1 2HGP -lift ♦ 1 p 1 4 

.24/6-12HGP-HG+1P1 5 

. 325F- 1 8HGP-HG4 1 p 1 6 

• 3 726 - l 2H5P-MG+ 1 Pi 7 

. I ft2 L - ] 3M6P-HG ♦ 1 P 1 1 

• 266F-1 3M(iP-MG+ 1 PI 2 

. 382F- I 3h(>W - HO ♦ 1 P 1 3 

. 5 3 1 6~ - 1 3HftP-HG ♦ 1 P 1 4 

. 7826 - 1 .UKjft-HG* 1P1 5 

. 9 1 2P- 1 3nG k-MG ♦ 1 P 1 6 

• 1 2 1 F - 1 2 • 1 ft 9 - M G ♦ 1 P 1 7 

.20 OF - 1 3M»*-HG‘*i ► 5 1 

, 22 / F - 1 JhG ♦ -H(jM ♦ 5 2 

• 25 U 6 - 1 3MG ♦ -MGW ♦ ft 3 

. 1 6 6*'- 1 3M(1 5 4 

• 9ft / 6‘ - 1 4hft ♦ -tiGM ♦ 5 b 

,6 1 .3F- 1 b 6 

• 4 1 ftp - | 'iMi ♦ -hGm ♦ 5 7 

. I 1 *.6 - 1 3H(, 4 -HGM# b 1 

• l 3UF -1 3H(* ♦-HG-m# b 2 

. 1 4 0 6 — J jhh ♦ - H ( i M 4 ft 3 
■ 14 6 6-1 .3 M ft 4 - HG M ♦ b 4 
. I 486 - 1 3 m ( 1 4 - H 1 , M 4 S 5 
. 1 4 7 F- — 1 JhO ♦ - HfjM ♦ b 6 
. 1 4 1 F - 1 *Ml. ♦ - HGM 4 b 7 
. 1 546 -1 iHf, ♦ -M(,M ♦ 3 1 

• 1-516—1 1 H ft ♦ - M ( > M 4 3 P 

• 2 M 06 -1 ^«ft 4 -r i(,3 4 3 3 

• 1 P)r - 1 1 lift ♦ - M ft • 1 ♦ 3 4 

• 7 »9F - 1 ♦ -HGM4 3 s 

• - 1 4 M(i ♦ - HftM ♦ 3 ft 

• 1 ’ft* 1 - 1 i f 

,6 8ft, - 1 '♦ Mft 4 - HGM 4 i J 

.8 1*6 — I 4 M « 1 » — M 1 1 M ♦ 3 2 

.9 17, - 1 4Mft ♦ - 4 3 3 

• i 0 1 1 - I Ini 1 ♦ -M6" 4 3 4 

• 1 (iftK - 1 Gift 4 - HftM ♦ 3 ft 

.10 76-1 INI, 4 -MOM , 3 ft 

• 1 856 -1 1M(. ♦ -H6G ♦ 1 7 


APPENDIX U 


The computer program "PROP" is listed below. It can be used 
to determine the values of the volume averaged plasma properties and 
the uniformity factors needed by the computer program "HG." Th° data 
needed to determine these quantities is obtained from a Langmuir probe 
survey of the discharge chamber in which the plasma properties are 
determined at many different locations within the chamber. This data 
is used to numerically evaluate Equations (10) to (15) and (21) yield- 
ing the volume averaged plasma properties and the uniformity factors. 
Comment cards are included in the computer program to indicate the 
purpose of each section. A CDC 6400 computer will use approximately 
thirty seconds of Central Processor time to evaluate five sets of data 
obtained from five Langmuir probe surveys. 
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program prop < t nrut * qu tpu t * r apfs = input * t a pe: 6 = output i 

This PROGRAM CALCULATES TmE AVERAGE pRDRFRTlh'S 
PL*5-PA0IAL POSITION OF LANGMUIR PROhF POINTS ( 1 -CF NTEHL INF ) 
PCS-POS. OK OFSTPKO OaTA POINTS IF NR # Nb# 4 

P0S7-AXT A L POSITION t)E LANGMUIR RROhF POSITIONS < 1-UPSTREAM PNT, 
NT-NO. OK TRACES Pb> SKT NR-NO. «)K RADIAL POINTS 

I F L AG- 1 If ONE WANTS To PRODUCE! A SET OF POINTS UPSTREAM 
IKLAG-? POINTS PRODUCED AT The RAFFLE 

DIS-DISTaNCF from SCREEN To POINT WHERE ThK GENERATED 5FT*IS TO 
PF PIACKO 

PCATH-CATHOOF RADIUS 

CNP-LEMGTH OK PRIMARY KLFCTNQN REMON AT CF NTtRL INK 

DIMENSION XNP(7?)« yr.:w(/P)t T<70>* 7</n j, v(soj* Util), u(21)t 
121)« TF(?1 ) • AMS) < RMS) • CCISJ* 0(3(5)* EE ( 5 ) * F F’ t S ) « FND(ll)* 
2S(S)» W(tl), ROS?(l))* YT ( tl i * Y P ( 1 1 ) • Ym(11>, y/(11), POP ( 10) ♦ 

3 A ( 1 0 ) * VF(S1) 

I N* F S T 0 M AT (70) * A / ( 7 0 ) * AN? (70)* ANM(70)* A(ll)* H(M)* C(ll)i 

KID* E ( 111 * VT(ll). PI. 5 ( 5 ) 

OATA AA*RR*CC*r)Q*FE. / 25*0 # / 

DATA POS /0. * 1 . 1 *?.?* 1. 3*4./ 

OATA PL 5 / 0 # * 1 

TNTFGRATFO CROSS SECTIONS FOR ♦ TO ► 

PFAD <5. 1?9) < 7 I ( I) *F ( T ) * 1 = 1 *21 ) 

WFAO ( • 1 P9) ( TF (T) * G ( I ) *1 = 1*21) 

WRITE (5*133) <71 (I) »F ( I) *1-1*21) 

WRI TK (5, 133) ( TK ( I ) t G ( I ) *1 = 1*21) 

pi =6.?83? 

101 REAP (5*t?<S> N R * N T * 1 F L Mm 0 I S * R C A T H * E NP 
IK (EOF(S)) 1?5*10?*125 
10? ND-NT/NR 
N N = N n 

IK < IFLAO.NF ,0) MN=ND*1 
READ (5*127) (POS? ( I > * f-1 *NQ) 

RADIAL DISTANCE TO CkITICAL FIELD LINK 


PRP 
PPM 
PRP 
PRP 
PRP 
) PRP 
PRP 
PRP 
PRP 
PRP 
PRP 
PRP 
per 
PR P 
F (PRP 
PORRP 
A A PPP 
PRP 
DPRP 
PRP 
PRP 
PRP 
PRP 
PRP 
PRP 
PRP 
PRP 
PPP 
PRP 
PRP 
PRP 
PRP 
PRP 
PRP 
PRP 
PPP 
PRP 
PRP 
PRP 


RKAD (5,127) ( F NO ( I ) * I - 1 * NN ) 

READ IN ThE PROPERTIES 
DO 103 I - 1 • M T 

103 READ (5*130) T ( I ) * 7. ( I ) * xmp ( I ) » XNM ( I) 

IF (NP.fO.4) tjO TO 10H 

I S - 0 

DO 105 1-1 *nT*NR 
NND= I ♦NR- 1 
DO 104 j=T.NNO 
K - .J ♦ 1 - I 
YT <«) *?T { J) 

YM (K ) * X NP ( ,J ) 

YM ( K ) - * N M ( J) 

104 Y 7 (*)-/( U ) 

DO lr>S 

l = . ! ♦ I S 

tail a t t * k n ( r»j , t t * n r * i * pc s ( j i * a r ( l ) ) 

TAII A 1 Ir EN (P| - * Yi> , 1 ,HO^ ( j ) , AIR ( I. > ) 

r A! ( A Mrf N C Ml -1 * r * ♦ Km , 1 « Kiis ( . t ) * ,\.V4 (( ) ) 

1 OS CALI AITKFN (P|„ StY?|-4M*1 ***()S ( j) . A/M ) ) 

10 5 T c - T ^ ♦ 4 

NT - \ T °4 / * R 
M'::NT/4 


PRP 

PRP 

PRP 

PRP 

PRP 

PRP 

PRP 

PRP 

PRP 

PRP 

PRP 

PRP 

PRP 

PRp 

PHP 

INP 

PRP 

PWP 

PRP 

PRP 

PRP 

PwP 

Pm P 
R t J P 
Pi'P 
P»*vJ 


10 
20 
30 
40 
SO 
50 
70 
80 
90 
100 
1 10 
120 
130 
140 
ISO 
150 
170 
180 
190 
200 
230 
210 
2?0 
240 
250 
260 
270 
280 
290 
300 
310 
320 
310 
340 
350 
350 
3/0 
380 
390 
400 
410 
4 ? 0 
4 30 
440 
45G 
4*0 
470 

4 R 0 
4R0 
500 
510 
520 

5 30 
540 
550 
550 
5 70 
5 M 0 
590 
5 0 0 
M0 
*20 
5 U> 
5 4 0 
5 SO 


r> n n o o r> 
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DC 107 J-UNT 

H 9P 

66 0 


r ( j ) . a r ( j ) 

PPP 

6 / o 


7 ( J) -A/ ( J) 

PPP 

6 ft 0 


XNP ( J) ~ANP (.J) 

PPP 

69 0 

107 

X^M (j) 

PPP 

/no 


VfftITF U * 13 0) ( f ( N ) • / ( Ti } » XMP ( N ) • XN M ( N) f N- 1 ♦ NT ) 

PPP 

no 

10ft 

IF ( IFlAG.F 0,0 ) CO TO 1 1 /« 

PPP 

/80 



pPP 

/J0 


CALCULATION op ThK pATPA Sf T OF POINTS 

PPP 

/A0 



PPP 

/SO 


N=1 

PPP 

760 


oo no r = if4 

PPP 

770 


L - 0 

PPP 

780 


no 1 no j=t *nt ,4 

PPP 

790 


l=i n 

PPP 

8 00 


Y T ( u ) a r ( J ) 

PPP 

ft 1 0 


YP (lj =KNP{ J) 

PPP 

880 


Ym (D sxNM ( J) 

PPP 

830 

109 

YZCL)aZ(J) 

PPP 

ft 4 0 


k =NT ♦ T 

ppu 

ft SO 


CALL AIT<FN (POSZ»Yr*Mn«N*OIStAT<K)> 

PPP 

ft 6 0 


K = NT* r 

PPP 

ft9Q 


CALL AITkfn (POSZf rTtM)«NtDIS*AT(K) ) 

PhP 

86 0 


CALL AITkFN (POS/*Y/«M)tNtDIS*AZ(K) ) 

P*P 

8/0 


CALL A I TkFN (POSZ t YP >fif)»NtOIS f AN-P (K ) ) 

PPP 

ftftO 


CALL AIT>KN ( POSZ • TP tNU* N » 0 IS* ANM (k ) ) 

PPP 

690 


IF ( AT (K) ,LT. 0. ) AT{K)=o.O 

PPP 

900 


tF ( A7 (K) ,LT • 0. ) \/ (* ) so ,0 

PPP 

910 


TF ( ANP <* ) # lT ,0. ) ANP(*)*U»0 

PPP 

980 


TF IfLM(K) .LT.O.) AhM(" ’ =0.U 

PPP 

9 3 0 


P I T F (S,1ZH) A T ( ^ )•«/{), ANP <k >, ANM <k ) 

r> K p 

9^0 

no 

CON T INUP 

PPP 

99 0 



PPP 

96 0 


PFShuFM TNG OF POINTS 

PP p 

9 / (j 



ppp 

9n IJ 


N 0 = N n ♦ 1 

PPP 

990 


nnn=no-i 

PPP 

1 0(10 


XK:P\D(M^] 

PwP 

n i o 


nc 111 J ~ , NMN 

PPP 

10 80 


J-NHO “I 

P P P 

10 80 


FMKJ) - b NO( J-l ) 

PPP 

1 0 '* 0 

1 1 1 

POSZ ( J ) - k’O 97 (J-l) 

PPP 

1 090 


F N D ( t ) = > « 

P P P 

10 6 0 


PCS7 ( 1) “T* I S 

PPP 

10 70 


NT = \T *4 

PPP 

1080 


NH-M/4 

PPP 

10 9 0 


DC 112 TsS.Nf 

PPP 

1 100 


K s N T ♦ S - I 

PPP 

1110 


J = K - A 

PPP 

1 180 


T ( K ) = T ( J ) 

PPP 

113 0 


?(*)-?{ J) 

P *P 

11^0 


XNP ( K ) = X6P ( J) 

P (J 

1 iso 

n? 

XNM {* ) = X**m ( J) 

PmP 

1160 


MmtnT-3 

11 

1W0 


no m I-NTm , NT 

Mm m 

1 lHO 



Pi- m 

1 I'M 


T ( J ) rr a f ( I ) 

( ►- P 

1/ Ml 


7 ( j ) - a ; ( i ) 

Mi. i. 

K 1 *i 


t P / , i > - A f.r* ( I ) 

* M ^ 1 

1 

1 1 j 

UM(.I) “A ru( { ) 

i ... ( 

1. 


ST Am t i j ( *u '• A( pm. a T T on' 


1 4 .. 

1 ' • 

i l '♦ 

no 1 1 t n .no 

? i, n 

Jit 


; -**(1-11*1 
* - 1 ♦ * 

* ,5 



His PHnmiM S AN A Hu AY OF KdM Af CONSTANT / HO I NTS 


PHP 


P»P 

C P HP 

00 11/ J = t. *K HHP 

m- j-t ♦ 1 PPP 

W r 1 3 0 S ( M ) « H I P H P 

CAM AT I h K N <1* f ( J) fl«I) P HP 

CA||. ATTKKN ( / I 1 WW I J) f K t ) P HP 

X T - XNH ( J ) ♦ KM-' ( J ) .'HP 

TF ( H m { j) »(■ ( J , 0 « 0 ) on TO 1 IS PHP 

TH.M r 1 , (j ♦ XMp ( j } /XNV ( J) PHP 

GO TO 116 PHP 

i is n hm,-. 1 , o php 

ns VF ( j) = r WT ( T ( J ) *4 .HU I tPM) PHP 

A A ( V ) - Y.\P ( J ) «XT»K I»H PHP 

HP ( **) -* X 1 «P PHP 

fT(W) =X[ »H*XNH (J) PHP 

pp ( m ) = jca.v ( j) <*X r <>0 1 *W PHP 

KF(M)=yNM(j)-»xT*M P IT P 

117 KF(MJ-H PHP 

C PHP 

c IMFOHATTON OK K m irMNCi G ( / > PHP 

C PH H 

PHP 

NM=lN PHP 

CALI T \' T *< ( HOS . A \ • ■* » '-4» i H * FNTH I ) #H) PHP 

A ( I ) = H p HP 

IK ( IM.AG.M’ • 1 .op. 1 .NF • 1 ) GO TO lib PhP 

r a M I \ r.. ( P0S« A A « '* fp * HCA TMtMM) PHP 

A ( I > * A ( I) -n" PHP 

11M CALI T nTw ( Pf’lS tW' ♦ - J * NP • ? « F NO ( I ) «H) PHP 

H(I)=h PHP 

IK ( I PI AC-. . NT . 1 . (».* , I . NK . I ) hlj TO l 1 W PH lJ 

TAM t \ I »■ { LOS ♦ M< * • NM , ,J , HC A Tm.mM ) PHP 

n ( I ) - p { I ) -i-M php 

110 CALL IS Tf { POS * ( . f • •' « V ♦ ? % K NO ( T ) » H ) P^P 

C ( I) =H PHP 

IF (IM A '■» , NF . 1 . OH * I . W . 1 ) HO TO 1 HO PHP 

CAM. TSTH (uOS»CC»‘-’Pi M i 3 »^CATHiHM) PHP 

C ( I ) “C { I ) -pm PHP 

1?0 FAIL TMiJ l KOs.fiM ’MNM NO { l ) » H) Pw P 

!)(1)"H PHP 

I F ( I K l A b . F F • 1 • f • 1 • OF ,] J GO T(j 1 ? 1 PHP 

CALL IN Tp ( lJ OS i Li; • '* • v* « ^ i kC A I M »HM) PHP 

n ( I ) -r* ( i » -*t" php 

1 p 1 CAM. TNTP ( Pi'S .MM '> , rgM , f J , | ; *j(i ( I ) *m) PHP 

F ( 1 ) :: H P U P 

IK ( T F i *. i o . M • 1 • i ■ ^ . I • • J 1 * 1 ) (iO Tm IrV PHP 

r A I l I S T P ( , > n s , f F • ‘ * t N • i* C A T h ♦ 1 1 i ) ppp 

F ( M - F ( I) -*<M PHP 

lam fsrp < ims ,pf ( m . *4*.*, / t kNi; ( I ) *m j php 

v r ( i j -m php 

IF (If! ‘ . \-F • 1 I • V . 1 ) no TO 1,M l 'HP 

CAM. 1ST* F K *v, v • -* t ' ; C A I • i*-’ 1 ) t; Hp 

v r ( i ) .. v r ( T ) “ J F- P 

r < \ t I M o p w ' li 


T Tv T 1 ' • A T ; . .v Li / (' I v M r [ (.% 


MOO 
lUO 
1 VO 

i no 

1 340 
i no 
l no 
l no 
l no 
i .no 

14 00 

14 10 
UPO 
UJO 
144 0 
14S0 
1 4 b 0 
U TO 
1 4 M 0 
1 4 HO 
1S0G 
1 S l 0 
ls?o 

15 30 
1 S 4 0 
1SS0 
hhO 
\^l 0 
1SH0 
lsou 

1 b 0 0 
IS1 0 
1 0 
1 nio 

1 s 4 Q 
1 SSO 
1 S b 0 
1 b70 
1 b8 0 
1 bHO 
WOO 
1710 

\r?o 

1 7 HO 
W4 0 
1 7 S 0 
1 7 b 0 
1 7 7 U 
WPO 
1 /GO 
1 m 0 0 
l MO 
1 VO 
1 'MO 
1 m/. t; 

1 - M 
1 > f i 
!■»/(, 
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CALI AITM'H (POSZ t V •hP*/ , *0 .*)? HP) 

OMOH-PH 

CALI OTfr, ( pos / * C « V t M) ) 

CALI AITKFN (PnS/ # VfNlH* 

CALL AITKKM (POSZ*V*NO** 

oc=qc-pc 

CALL OTFG (POS7#Hf V«ND) 

CALL ATTkFN ( P(*S Z * V • NU • * 

CALL AITk-FN {PnSZ.VtNOt* 

on=on-Pn 

CALI. OTKG (POS7iP'«V*Mn) 

CALL AI T*FN (POS/*V*M)t* 

CALL AITKFN (PPSZ« V *MUt * 

UF=GF- PL 

CALL QTFG ( POSZ * V T « V « NO ) 

CALL AITkfm (POS7*V*NM*? 

CALL AITkfN (PPSZ*V*Mi.M* 

QVT=OVT-PVT 


?*FNP*QC) 
! * 0 # 0 • PC ) 


• tNP*urj) 

•O.Ot«D> 


fFMPfOF) 
. O.OiPF) 


•EMP«uvr> 

«0.0»WVT> 


1*4 

1?S 

IPS 
1*7 
1*B 
1 *9 
1 30 
131 

1 3 * 


CALCULATION nK PPOPf>TltS 

F Z I =0 A /OC 

isTE*DO/OF 

A\FsO«/OVT 

AP=CC/OP 

A^-OK/On 

0=QC/0F 

x p = sop T ( ( ap <> r i ) / ( 1 . ♦ l • /u) ) 

XP=*OWT { ( A>" Af.f ) / ( 1 . ♦ i; ) ) 

XNF=XP*XM 

CALL A I TKP K ( * IF «*) . 1 .0 rt » tFS) 

C A I l AITktkj ( F . Z t * * 1 » 1 ♦ F / I » Z I S ) 
w i x P / X m 

a P I T F (MllO OA .op . oc * on t UH 

w^itf (p.nn * ^ * x « r *■: s * 7 1 s 

^ Pc I T K (F . 1 V) XNP . 0 V r 

Vf- A = S(JPT ( TFPi»/* l f«fUF ( 1 . *P) ) 

F I NO NOMI IN I F t)*/M I T Y K AC T OPS 

CALL APh A { XF»p * XNM • XNt • H(.»5 %POSZ »NP t *<LATM « VF « Vt A ) 

IF ( IFl.AG.FO, 0) GO T() 1*4 

NT-NT-4 

NO-M/4 

CONTINUE 

GO TCI 101 

STOP 


f-OP^AT ( 3 T S . 4 F 1 0 . 4 ) 

FOOmaT (FPW'.S) 

FOPM aT ( 1 0> * '»*■' 1 0,4 « *F 1 0 . 3) 

KCPpaT ( 3 ( V 1 (1 • * * K 1 0 , 3 ) * * 0 X ) 

fowmaT ( l ox »?m o ,4 i n , 3 ) 

LCPpaT <*X« 1 l MNP ICM-U Sifc!l,<M lx, lO'NM ICm- 3) =*Fll t 4«JX* 

1 P/N*^“ , F 1 o .4 * 1 A f H.Ht (r V) -«F1C.4*3X* F N . ( F V ) = . F 1 0 . 4 ) 

row at nox. p‘*m k.t, i *m l *4 tsx • icm n r , r 1 o .*♦/) 

FOHM at < /* (r u .4 t *X ) ) 

F NO 

MOPOLTpF P-Tc. ( HP *M\«M ^>-'0^ I'jTiXf'pl') 

T^IS MU’-.'Ht tGI IMm'-UIS Y(X) AML' NF'< FINOS ThF V A 1 l j F OF Tm|- 
t \ T K C- P A | AT X >■ ■> 

PT^N-STCN it ( T m ) • YtM'ii i A ( 1 to • u ( ) •' M < r < 100) 


ppp 

10G0 

ppp 

19 7 0 

ppp 

19P0 

PPP 

t 90 0 

PPP 

*000 

PPP 

P010 

PPP 

*0*0 

PPP 

*0 JO 

PPP 

*040 

PPP 

*oso 

PPP 

*oso 

PPP 

*070 

PPP 

*0rt 0 

PfJO 

*090 

PkP 

*100 

PPP 

*110 

PPP 

*1*0 

PPP 

*1 JO 

PkP 

*140 

PPP 

also 

PPP 

*1G0 

PPP 

*1/0 

PPP 

*1«0 

PPP 

* 190 

PPP 

**00 

PPP 

**10 

P^P 

***G 

PPP 

**30 

PPP 

? * 4 0 

PPP 

?*so 

f > P p 

* * G 0 

ppp 

**70 

PPP 

?*M0 

H r* P 

**9U 

P v p 

* »oo 

PPP 

* il 0 

PPP 

* 1?Q 

PPP 

* 3 3 0 

P PP 

*340 

PPP 

* JG 0 

PPP 

* 3 G U 

PPP 

* 370 

PPP 

* 3B 0 

PPP 

* 390 

ppp 

*400 

PPP 

*4 10 

PPP 

*<4*0 

PPP 

*4 30 

PPP 

*440 

PPP 

*h-»0 

1- i* P 

* '♦ G 0 

i • UP 

*4 7 0 

I'tip 

* 4 P 0 

PHNF pp 

* 4 <i 0 

I'f. p 

*SO 0 

PI. p 

*S 1 0 

p.cri 

* G * 0 

pi.P 

*G 3 J 

1 l«l> 

* G *♦ G 

f. T , ; 

1 L 

M Mi 
f f r. 

*0 

MM. 

4 ■ J 

M r t, 

G 'J 

r 

* 



--1 


A 


' 


-> r\ 
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CALL MMH'P ( » * Y • A . N I F tN'JSI- ■ I') t IN i> N10 

CM.I. H T m lAiHfC^'ir-ni) NTS 

CAM A M k F N { A f r \j t l ► r • Xp I YP ) NTS 

PF TliPN NTS 

c NTf » 

KND NTS 


SUOPOUT INF APF A (Xr?M * XM|.| t XNLF f')TP*rTYtNLNt p » Vht V£ A) 


C 

r 

r 

c 

c 

c 

c 

r 

c 

r 


fHl< MJppHJ T p/f CAiruiATFS fHF Pl.ASMA 
PPY-Ptilf-TS (/) K IMi PKlMAPV KfM.O 

PPP-KOIMS (*') Nt K I N I NG P ri I »> A p X HM.H 

F; P P F -NO • nr p.<I*MWY F I ► L ( l INF POINTS 
I P 1 T F L nSSF S APF ALL (t'jfn 1 0 ThF 

Cl* 0 I T A N f F Tn 4 A F F i F FP<JN SCPKFN 
C7U-POIM WhKPF PPlKAPY f IKU) LINK iMTFPSKCTs 
(IF I FL A G = 1 ) 


fiONUNlF 0«?MI fy FACTORS 
l I Nh. ( F P*AM Sf'PFFN OK 1 1) ) 
LINK (FpOM OnTF.w LINK) 


CATmjbt POLK PIFCF. 

CATHOOK POLK PIFCE 


o I MF'NS 1 OM n TP ( * ) t f»T Y ( 1 0 ) 
1(10)* P^Y ()0)* p ( S ) * / ( “> ) 

?K 1 ( 1 1 J t VF P ( 11 ) « v» Jlllli 
PF AL II * r ? 

N F I.. A G - 0 

M F L A G = 0 

h ah ( s , i o n 

C^o-n I 

pk A n (Still) 


XT 1 ( 1 1 ) 
ANP (Si) 
VK4 (11) 


x \? ' n > 

XNM (SI ) 


M'KlDt X I 4 ( 1 1 ) * 


APA 
AKA 
AW A 
AKA 
A PA 
APA 
APA 
APA 
APA 
APA 
APA 
PPP APA 


* 0<S> t VK (Si ) ♦ V (S) t 


NPPF tin A G • L I *CPU 

( K P P ( I ) * P P T ( I ) » I - l * NPPK ) 


Mf[, t \ ( A I Cl'l AT TON'S 
SHjFF l F K'n „ p. F MS I t 1 K S 


into coppkct a prays 


on iul lit op * 4 
J=M 
< ~ T ♦ 7 
\ I ♦ '1 
M t I / A ♦ 1 
V F 1 ( M ) •' V F ( I ) 

VF r ( “ ) “VF ( ,J ) 

V * 1 ( ' ‘ ) - V F ( * ) 

VK 4 ( ►<) -VF' (|. } 

* T P (f‘ ) - AMM J> ♦XNM ( J ) 

X T 1 ( ♦• ) - XNR ( k ) ♦ AfjM { « ) 

X I 4 ( *M P ( L ) ♦ XL-!-* (L ) 

101 X 1 1 (M ) t \.p { I ) ♦ XNM ( I ) 

*KITF OhT TmFSF APPAYS 


W 1 TP { 

* P 1 T F ( 

^ - T r r ( 

f. k i r k ( 

a / n . 

i Y - <■' r * Y ( 

A L - «) . 

a k - n , 

Y / ? 


. 1 1 0 ) (MI(T)O--l • NL N ) 
. 11 (i > ( X I? ( I ) * l = 1 tKLN) 
.110) ( X I M I) * f - 1 • K'l N) 
.110) ( X 1 * * ( T ) • I = 1 tNLN) 

p F ) / s 0 f ) , 


Y/ 


* I A i f I S 1 F ‘ii’ A | 


- I “ A 


f N 1 


• I 1 -Pi 


1 • Ml 1 


VAPA 
APA 
AKA 
APA 
AKA 
APA 
APA 
AKA 
APA 
APA 
APA 
AKA 
AK A 
APA 
AKA 
APA 
APA 
AKA 
AKA 
APA 
APA 
Ak A 
APA 
A A 
APA 
AKA 
APA 

• A 

• A 
PA 

i r» A 
ip A 

• A 
AKA 
A >« A 
A i « ^ 
Am A 
Am A 

. A 


A I* 
A k 


AK 


A 

A ■ « 
A (■ 

A • 
A ' 
A ■ 
A ' 


K0 

no 

100 

no 

120 

10 
P0 
30 
4 0 
SO 
SO 
7 0 
K0 

no 
100 
110 
1P0 
1 30 
140 
ISO 
ISO 
170 
IhO 
l no 
POO 
P10 
PPO 
P 50 

P 4 0 

psu 

pso 

P MJ 

f*:;0 

P)0 
3 0 0 
3 1 0 
VO 
i )U 
14 U 
ISO 
JSO 

3 ro 

.3K 0 
3SV 
'♦0 0 
’♦10 

4 f 1 0 
4 10 
44 0 
4 K 0 
4 h 0 

4 7 0 

4 M (' 

•IPO 

■> 1 (. 
•) 




t / s 


o r> o 


CALL AflKFN ( PP Y f PuP • FPPK f l , T p , |)p ) 

('ALL MT*f\N (PPY t PPP tNPPF f 1 « YM§OM> 

ox = np-r)M 

N-l 

IE <T.GT # ??5) N=1 

CALL A I T K K N (DTY* X I 1 »MLN*N* Y£ »o ( 1 ) ) 

CALL AITKFN (nTY«XI?*M.MiNiYZ # P(i?)) 

CALL A 1 Tk FN (D r Y » X n f MN t Nf Y/ ,P (3) ) 

CALL AlThFN (OTY*/I<*iNLNtNtYZ.P(4)) 

CALL C*FCK ( P i 4 ♦ N F L A fj ) 

J = 4 

P(5)=P(4) /5, 

CALL A I TXFN (OT Y « VE 1 iNLN*N* YZ t V ( 1 ) ) 

CALL ATThPn <nTY«VF?,NLM*N«YZ*V{Zi) 

CALL AITkFN (GTY « Vf. 3*A*LM*N« YZ« V (3) ) 

CALL AITkFN { 0 T Y * V F. 4 * F L N * N t Y Z t V ( 4 ) ) 

CALL ChFCK (V*4*MFLA0) 

V(5?=V(4) /8. 

CALL AITn-FN (OTP.P, j,k, xp, XNE) 

CALL AITkFN ( OTP ♦ V t J • N * XP * VEE ) 

YZ = Y Z ♦HY 

AP=AR4S0PT(nx**?*0Y**?)*6.2H32*XP 
A/=A74SO*T (DX**2W1Y<>*?) *h,?R3?<*xP*XNF <►*/?* VF.F 
10? A N = 4 N ♦ S (; P T (nx««?*nY**2)**.?832<*XP*XNE 0 VhE 
C 

C SECTION fop GPIOS 

c 

p (l)=*>.?P1? ft 0TP(l)«XIl(NLN) 

P(?)--6.?R3PonTD(?)#XlP(NLN) 

■=> <3> =6.?«3?*r)TR (3) *X I 3 (N( N> 

p t 4 J s6«?Ml?»f)TP(4)«Xl4(Nl N) 

Call chfck (p*4*mflag) 

P (S) sO, 

V ( 1 ) =VF ] <NL^) 

V ( 2 ) = VK2 ( N I N ) 

V ( 3 ) =VF3 (Nl N) 

V(4)rtVE4(NLN) 

V (*> ) sQ # 

DC ini ik=i*s 
10 3 0(lK')-V(lK)*P(fK) 

CALI OTFCi (0fP»0*Z*b) 
wPITF (A*lOh| 7(5) 
n = ? (Sj 

A N = A N * 7 (5 ) 

00 104 Jk=2*5 

104 Q(I*)=Q(TK)*P(TK)/(f *?h3?*0TR(JK) ) 

CALL OTFCi (DTP*Qt/«S> 

1 2 = 7 ( 5 ) 

A7=A?*7 (P) 

7 (5) =7 <S) 

*PITF (8.108) 7(5) 

AR-=AQO.l4l59*OTP(i)**7 

0 (5 J =0. 

SECTION FOR PAFFl h 


N - 1 

CAU a I r*FN { f) T Y t X I 1 ♦ M N • N « 0 I « P ( 1 ) ) 
CALL A I T k C f j (oIYitt^M f|«N.OI»F(F)) 
CALL A J T * ► N (tUY«xp»MfuM.ni»tMi)) 

Cau. attkfn { o r y « a i >* * *[ n « n * o i # 4 j ) 

CALI. CnfC^ ( w f 4 « ** L A 11 ) 

CAL! ( 0 T Y • Vf 1 * M N * N « 0 I • V ( 1 ) ) 

CALL A I T « f N ( O T Y « vt e • \| N « N • 0 I « V ( 2 ) ) 
CALL A I f k F N ( ! * I Y « v M • M \ « W • U I • V ( 3 / ) 
C A l | AfT*FN |f'TY,vf 0.1 '.i‘m. 11 «v (4) ) 


A PA 

*00 

APA 

*1 0 

APA 

*20 

APA 

*3 0 

APA 

*4 0 

APA 

650 

APA 

*60 

ARA 

670 

APA 

680 

ARA 

6 9 0 

APA 

zoo 

APA 

no 

APA 

7r?0 

ARA 

730 

APA 

7 4 0 

ARA 

750 

APA 

76 0 

APA 

7 70 

APA 

7 8 0 

ARA 

790 

APA 

600 

APA 

810 

APA 

820 

APA 

630 

APA 

840 

APA 

850 

ARA 

8*0 

APA 

370 

APA 

880 

APA 

690 

APA 

900 

APA 

910 

APA 

920 

ARA 

9.3 J 

ARA 

94 0 

APA 

980 

ARA 

9*0 

APA 

970 

APA 

980 

APA 

9 9 0 


ARA 1000 
A PA 1010 
A I? A 10 20 
APA 1030 
APA 1040 
APA 1050 
APA 10*0 
APA 1070 
A P A 1 U R 0 
APA i ova 
APA 1100 
APA 1110 
APA li?0 
APA 1 1 10 
ARA 114 0 
APA 11S0 
APA 11*0 
APA 1170 
A*a 11 PO 
APA 1 1 NO 
A p A ? > 0 0 
APA 1 2 1 (J 
APA 1,VU 
APA 1 2 3 0 
A p A l ? 0 
APA 1 /SO 


o r> o o o o 


80 


CALL CHECK ( V.4 »MF| AG) 

0 ( 1 ) = fi.?A‘*?*OTR ( 1 > <1 ) *V < 1 ) 

0(?)s6.?P12*OTP(2)*P(2)»V<?) 
Q(3)=ft.76 37'MUP(3)<>P(3)*V<3) 

0 (4) =6.7832*0 TP (4 ) *P (4) *V ( 4 ) 

CALL OTFG (PTrff')t7tS) 

CALL AITkFN (HT«t St l.rtttC) 

DC ins I J = 1 ♦ 4 

105 Q(IJ)=Q(TJ)*P< JJ) 

CALL OTEG (0TPtQ*/*6) 

CALL AITkFN COTP, ; *S» 1 ♦*•£0) 

AP = A«0.141«5****2 

AN=AN*EC 

A7=A7*ED 

IF ( IFLAG.NE. 1 ) GO ro 1 0f% 

CATHODE pole piece section 

CALL ATTKEN (DT r I 1 .ML N*N,C7UfP ( 1 ) ) 
CALL AIT*FN (OTYtXl?tNLN*N.CZUtF<2) ) 
CALL AITKFN < OT Y * X M * NLN • N . C 7U . P ( 3 ) ) 
CALL AITkEN (OtY : M4.\LN*N»C?U*P (4)) 
CALL C^KCk (P. 4, NFL. AG) 

P (5) rP (A) /H. 

CALL AITkFN ( 1>T » . P * - % N • P . X N I ) 

CALL AITKFN ( OT Y ♦ VE 1 ♦ NL N ♦ N . C ZU ♦ V ( 1) ) 
CALL AITkFN ( n T Y * VF ? . M.N . N . C 7U ♦ V ( 2 ) ) 
CALL A I TK Fm (HT Y . VM *M.N *N *CZU. V (3) ) 
CALL AITkFN < 0 T Y t Vt 4 ♦ M. N t N . C ZU * V ( 4 ) ) 
CALI. C HPfK (V«4«MFl.A0) 

V (5) =V <4) 75, 

CALL AITkFN ?PTw.V.5.N.P*VNl) 

CALL AITkFN < i j7 Y , X II • ML N ♦ N ♦ C ZO ♦ P ( 1 ) ) 
CALL AITKFN (UTV.XIE.NLN.N.CZD.P(^)) 
CALL A I 1 K F N (IJTY. X MtN|.N.N.C7t) *» J (3) ) 
CALL AITKFN (OTY«Xl4*Nt_NtN.CZD»P(4) ) 
CALL CHfCK (*>«4fNELAG) 

P (5) =p (A) /6. 

CALL AITkfn ( OT w . P t 6 « N ♦ P , XN2 ) 

CALL AITkfn ( OT Y ♦ VF 1 * NLN . N • C 71) * V ( 1) ) 
CALL AITkfn ( OT Y * VF 7 % Nl.N . N « C 7D % V ( 2 ) ) 
CALL AITKfN <uT Y . VI- ‘■UNLN.N.C7D. V (3) ) 
CALI AITkfn ( n T * ♦ vt 4 ♦ nl N • N ♦ C ZD . V ( 4 ) ) 
CALL ChFCk <v/,4,MFLAG) 

V (S) sV (4}/6. 

CALL AITkfn (DTR.V.5*n#p.VN2) 

XL-C/U-C70 

A P = A P ♦ X | *6.?83?«P 

XN= ( XN1.XN?) /7. 

VT = { VN1 ♦ V N 2 ) /7. 
AN‘iN*T,?^V»wM|^XN»V^ 

A 7 - A 7 ♦ * .7637<>h*L <> XN««7*VN 

wW I TF nijT pf GUI TS 

106 !f,im M 
v» f I T F (6.117) AN 

I T F (* t 1 1 4 ) A 7 
M=M>«»4i-/A‘j>VFA 

»F I 1 *XNEE ) 

W p I T f (^ • IP'1 M . *■ 7 « M w w 
wh I T F ( ^ ♦ 1 1 * : r,F l AG 

k p I T r ( f ♦ 1 f.P *, | AG 

w ► Tl.PN 


AW A 

1260 

APA 

12 7 0 

APA 

1280 

APA 

1290 

APA 

1 300 

APA 

1 310 

APA 

1 320 

APA 

1 330 

APA 

1340 

A»A 

1 350 

APA 

1 360 

APA 

1 370 

APA 

1380 

A P 

1390 

Ai- A 

1400 

APA 

1410 

APA 

1470 

APA 

1430 

APA 

1440 

APA 

1450 

ARA 

1460 

APA 

1470 

APA 

1480 

APA 

1490 

APA 

1500 

APA 

1510 

APA 

1520 

APA 

1530 

APA 

1540 

APA 

1550 

ARA 

1560 

APA 

1570 

APA 

1-380 

APA 

1590 

APA 

1600 

AHA 

1610 

APA 

1670 

APA 

1630 

APA 

1640 

APA 

1650 

AW A 

1660 

APA 

1 cj 7 0 

APA 

1 680 

APA 

1690 

APA 

1 700 

APA 

1 710 

ANA 

1 720 

APA 

1 730 

APA 

1 740 

APA 

1 750 

APA 

1 760 

AkA 

1 7 7 0 

APA 

1 / MO 

APA 

1 7 NO 

APA 

1 800 

AP A 

1 81 G 

A w A 

1 870 

APA 

1 me 

APA 

18*0 

APA 

1 85 c 

APA 

1 0 

A >< A 

1WC 

A p A 

1 68 C 

A wA 

1 8 p •. 

APA 

1 *o r 

APA 

lilt. 


81 


107 
1 OH 
109 

no 

in 

n? 

1 1 3 
114 

I IS 

II * 


F OK //. T 
FORMAT 

format 

1 t M I: S ) 

F f K m a T 
FOKmaT 

forma r 
f r w m a T 

FORMAT 

Format 

1.4) 
FCK^AT 
1 IMF? ) 


( ? I S * 4 F 10.4) 

<?*• VhI ir«.||:,Hl ( »| 

(?X, 4()HNM;aT IVt VM.OCITY EXTRAPOLATION OCCDWPEO . I 5 , 7 h 

(lOXfH (K 1 0 . 3 a ? A ) ) 

(HFIo.h) 

(37a 1 1HAN ( /C M .3 ) - a F 1 1 * 4 ) 

( 3 X a ‘*HAU ( ) ~ . t 1 1 ,4 ) 

(3Xa 1JHA/ (l>V/C^ ) - *F 1 1 .4) 

<?X. 3hF1 = ,Uii. 3HFa=,F 10.4.2X, 13HI»*/I* FACTOR »F 
(?Xt 3'-JHNFGATI VF TENSITY F X TPAPfJL AT ION OCCUPIED* 15* rtH 


ANA 
AKA 
TAKA 
AKA 
AKA 
AKA 
AKA 
AKA 
AKA 
1 0 A K A 
AKA 
Taka 


c 


AKA 


F:.r, 

St'-<K OUT I f'F CHh CK (K aN aNFI A(i) 

AKA 

AKA 


0 r MF NS I OH P (N) 

C K 


or. ioi i - 1 f ,\» 

CK 


T F f'MD.^.O.O) GO U) \ 0 1 

CK 


k ( r > -- o . o 

CK 


NFI.AfiafsjFL AG* l 

CK 

101 

FONT TMjF 

CK 


rftnpn 

CK 

c 


CK 


D 

CK 

« 


CK 


r^o 

1 '#30 
1940 
19S0 
i96 0 
19/0 
1980 
1 990 
?UQ0 
?010 
? 0?0 
? 0 30 
? 0 4 0 

?oso 
? 0 6 0 

10 

?0 

30 

40 

SO 

60 

70 

60 

90 

100 


